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PREFACE 


In  December  1975,  the  Automatic  Data  Processing  (ADP)  Center  of  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg,  Miss., 
submitted  a proposal  to  implement  and  evaluate  interval  arithmetic,  a software 
system  for  digital  computer  numerical  analysis,  on  the  Corps  of  Engineers' 
primary  engineering  computer — the  WES  Honeywell  G635.  The  proposal  was  later 
expanded  to  include  the  implementation  and  evaluation  of  an  interval  arithmetic 
software  package  on  six  different  computer  systems.  Engineering  and  scientific 
data  problems  were  selected  to  be  used  on  each  of  the  six  computers  with  the 
interval  arithmetic  software. 

The  work  was  funded  by  the  Office,  Chief  of  Engineers,  U.  S.  Army, 
through  the  Integrated  Software  Research  and  Development  (ISRAD)  Program, 

AT11,  Engineering  Software  Research. 

This  is  Report  5 of  a series  entitled  "Implementation  and  Evaluation  of 
Interval  Arithmetic  Software."  The  other  reports  to  be  published  in  the 
series  are: 

The  State  of  the  Interval:  Evaluation  and  Recommendations 
The  Honeywell  MULTICS  System 
The  Honeywell  G635  System 

The  IBM  370,  DEC  10,  and  DEC  PDP-11/70  Systems 

This  report  was  written  by  Mr.  David  A.  Cohn,  Mr.  J.  Brian  Potter,  and 
Dr.  Myron  Ginsberg  of  the  Department  of  Computer  Science,  Southern  Methodist 
University,  Dallas,  Tex.  Their  work  was  performed  under  Contract  No. 
DACA39-77-M-0105,  dated  24  February  1977,  and  under  National  Science  Foundation 
Grant  MCS  77-06611.  The  work  concerned  implementation  and  evaluation  of  an 


Report  1 
Report  2 
Report  3 
Report  4 


interval  arithmetic  software  system  on  the  CDC  CYBER  70  computer  system. 


Dr.  J.  Michael  Yohe,  Director  of  Academic  Computing  Services,  University 
of  Wisconsin-Eau  Claire,  developed  and  wrote  the  interval  arithmetic  software 


a 


j 


package  which  was  Implemented  on  each  of  the  six  computer  systems.  Dr.  Fred  D. 
Crary,  formerly  with  the  1).  S.  Arms  Mathematics  Research  Center,  University  of 
W1  aeons  In-Mad  Ison,  developed  and  wrote  the  AUCMF.NT  precompiler  which  was  im- 
plemented on  each  computer  system  as  a front-end  to  the  Interval  arithmetic 
software  package.  Dr.  Yohc  and  Dr.  Ctary  are  specially  thanked  and  recognised 
for  their  technical  contributions  and  assistance. 

Mr.  .lames  R.  Cheek,  .lr.  , formerly  with  the  ADI’  Center,  WES,  provided 
Initial  impetus  and  guidance  for  the  project.  Mr.  Fred  T.  Tracy,  ADI’  Center, 
WES,  provided  expert  advice  and  technical  guidance  during  the  project.  Dr.  S'. 
Radhakr Ishnan,  Special  fechnlcal  Assistant,  ADD  Center,  furnished  technical 
guidance  and  general  project  supetvlsion.  The  project  and  the  report  were 
monitored  by  Mr.  William  I..  Boyt  under  the  general  supervision  of  Mr.  D.  L. 
Neumann,  Chief  of  the  ADP  Center. 

Directors  et  WES  during  the  prelect  and  the  preparation  of  the  report 
were  COL  C..  H.  Hilt,  CE,  and  COL  .1.  L.  Cannon,  CE.  Technical  Director  was 
Mr.  F.  R.  Rrown. 

Copies  of  the  other  reports  of  the  series,  computer  listings  ot  the  in- 
terval program  and  ot  AUGMENT  tor  each  computer  system,  tnd  runs  of  the  bench- 
marks for  each  computer  system  may  he  obtained  from  the  ADP  Center,  WPS. 
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IMPLEMENTATION  OE’  THE  'AtdMENV'  PRKCOMP U 1.K  AN1>  INTERVAL 
ARITUME  Til'  OS  rill  COC  CYBER  70  system 


l.  Cyber  70  and  Oper.it  ln£  System  Arch  1 1 ec_t  ur  e 

Architectural  differences  amongst  the  6000  Series,  7600,  Cyber  70 
Series,  ami  Cyber  170  Series  are  sufficiently  minimal  for  this  report  to 
refer  unambiguous lv  to  CDC  computers. 

On  Control  Oat  a computers,  I/O  and  operating  svstem  functions  are 
resident  in  the  ten  to  twenty  peripheral  processors  (referred  to  as  PPl’a) 
rather  than  in  tin1  central  processor.  There  are  no  direct  comnninlcat ions 
between  the  CPU  and  the  FPUs.  The  CPU  cannot  be  Interrupted;  requests 
for  operat  ing  systems 1 and  1/0  functions  are  detected  bv  a PIT  polling  loop. 

All  communications  between  the  CPC  and  the  external  devices  are  through 
main  memory,  where  all  user  code  is  resident.  The  use  of  PPl's  for  1/0 
and  the  high  memory  bandwidth  are  kevs  to  the  high  system  performance. 

Extended  Core  Storage  (ECS),  if  present,  is  often  used  as  a fast  swapping 
drum  and  for  small  user  tiles.  Some  of  the  relationships  amongst  the 
architectural  components  are  illustrated  in  Figure  l. 

The  CDC  KRON'OS  Operat ing  Svstem  which  is  used  on  our  host  machine  Is 
highly  modular.  It  Is  based  on  the  Pile  Name  Table  (ENT).  Each  control 
card  that  Is  encountered  causes  the  ENT  to  be  searched  for  a user  file  to 
be  loaded  and  executed.  If  the  sequential  search  falls,  the  ENT  is  searched 
again  for  a PIT  program  or  system  CPC  utility;  if  the  second  search  suc- 
ceeds, that  system  routine  is  executed. 

The  Southern  Methodist  University  Computer  Center  has  a Cyber  70, 
model  77  machine  with  11  peripheral  processors,  no  ECS,  and  a maximum  main 
memory  user  space  of  approximately  200^K  words  (l .e. 6-tj^K  words) . The  Fortran 
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Figure  1:  CPC  bOOO  and  7000  Series  Systems  Configuration 


compiler  is  FTN  (version  s.b).  At  the  time  of  the  initial  Implementation 
of  AUGMENT  anil  INTERVAL,  the  KRONOS  2.1.0  version  of  the  operating  system 


was  in  use.  After  the  completion  of  the  Implementation,  the  Computer  Center 
changed  from  KRONOS  2.1.0  to  KRONOS  2.1.2  which  Is  compatible  with  NOS  1.0 
usej  at  other  CDC  Installations.  It  was  found  that  this  change  was  trans- 
parent to  AUGMENT/ INTERVAL , due  to  the  upward  compatibility  inherent  In  the 
CDC  product  (both  hardware  and  software).  Inasmuch  as  no  changes  had  to  he 
made  either  to  AUGMENT  or  INTERVAL  t ' make  use  of  CDC's  FTN  Compiler  or  to 
allow  for  either  version  of  the  operating  svstem,  it  is  our  belief  that  no 
changes  would  he  needed  If  AUGMENT/ INTFRVA1.  were  to  be  Implemented  at  another 
CPC.  site.  None  of  the  code  makes  anv  calls  to  specific  system  routines. 

Note,  however,  that  calls  to  certain  UNTVAC  Fortran  routines  (e.g.,  the  FLP 
function)  had  “o  he  replaced  bv  In-!'. no  code  to  achieve  the  same  purpose, 
and  rhat  certain  VNIVAC  routines  (e.g.  URW'ALK,  CRRT,  STRACF’'  had  to  he  re- 
placed bv  equivalent  routines,  coded  at  SMI'. 

We  did  find  that  AUGMENT,  even  when  overlaid,  was  too  large  to  be  rur. 

as  a tlne-sharlng  Job  (SMI  TELEX  permits  a maximum  Job  sire  of  lOO^k  words. 

i.e.  12.  ,K  words',  so  that  we  were  forced  to  run  In  hatch  but  this  was  a 
10 

minor  inconvenience,  at  worst.  It  is  our  belief  that  a host  computer  con- 
figured to  allow  larger  TELEX  Jobs  would  not  he  faced  with  this  problem  but 
it  is  not  necessarv  to  have  this  conf igurat ion. 
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1.  imp  1 ctren tut  Ion  ol_A!GMlNI 

AUGMENT  is  a Fortran-based  precompiler  that  facilitates  the  use  of 
non-standard  constructs  In  Fortran  programs.  It  can  also  be  used  to  aid 
in  translation  between  various  dialects  of  Fortran,  or  to  implement  non- 
standard operators.  It  accepts  programs  written  in  a user-defined  exten- 
sion of  Fortran  where  the  extension  consists  of  new  data  types,  new  opera- 
tors, and/or  r.ew  t.  library  functions.  It  requires  a description  of  the 
extension,  which  consists  essentially  of  a deser  ion  of  the  user-supplied 
routines  tl.at  Implement  the1  extension.  It  then  uses  this  description  to 
translate  the  extended  Fortran  into  standard  Fortran,  with  the  non-standard 
constructs  ttanslnred  into  references  to  the  routines  in  the  supporting 
packages . 


In  the  SMI'  Implementation  of  AUGMENT,  we  found  that  the  non-overlaid 
version  required  1 70oK  words  (i.e.  60,  K words'  of  core  memory  (and,  of 

o xU 

course,  a disk  channel).  Vsin’  the  overliving  suggestions  given  in  [O’, 
we  overlaid  AUGMENT,  which  reduced  the  size  requirement  to  llO^K  words  (i.e. 
40^K  words'  of  central  memory  but  at  the  cost  of  increasing  the  execution 
time  required.  We  were  able  to  follow  the  suggested  overlays  with  the 
exceptior  that  we  were  forced  to  move  one  routine  from  the  suggested  overlay 
to  the  root  overlay. 

On  CDC  computers  the  output  from  AUGMENT  can  be  compiled  on  the  FTN 
compiler  with  any  desired  compiler  options.  The  resulting  object  decks  are 
then  loaded  with  libraries  INTUB  and  BPA1.1B  which  contain,  respectively,  the 
routines  defining  Che  interval  and  BPA  operations,  and  the  primitives  which 
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perform  the  basic  arithmetic  operations  for  the  entire  package.  The  resu '.t 
of  the  load  can  then  either  be  executed  or  written  to  absolute  central  pro- 
cessor overlays  (see  Figure  2).  The  control  card  sequence  for  the  use  of 
the  AUCMKNT  precompiler  as  well  as  the  entire  Job  deck  structure  is  dis- 
cussed in  Section  6 of  this  report. 


AUGMENT  Source  Program 
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2:  Processing  Sequence  for  Program  Using 
Interval  Analysis  Package 
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1.  Implementation  ol  Imcwn  l Routines 

The  INTFRVAI  set  ot  routines  Implements  interval  arithmetic.  Interval 
arithmetic  is  an  arithmetic  system  In  which  rules  are  defined  to  create  in- 
terval versions  of  the  primitive  arithmetic  operations  and  then  to  combine 
these  intervals  to  create  interval  extensions  of  user-defined  mathematical 
expressions  and  functions.  The  application  of  Interval  analysis  to  a com- 
putational problem  results  In  the  generation  of  a rectangular  patal ' elpiped 
which  encapsulates  the  exact  solution.  In  our  implements t i on , all  the  usual 
arithmetic  operations  and  Fortran  library  routines  have  interval  counter- 
parts; in  addition,  several  Interval-oriented  operations  are  also  available. 
The  reader  unfamiliar  with  the  details  of  interval  analysis  should  read  the 
explanation  in  Appendix  R and/or  should  consult  some  of  the  references 
given  In  the  bibliography  in  Appendix  A. 

Since  chained  computations,  especially  those  Involving  multiplications 
and  divisions,  cause  the  size  of  the  resulting  intervals  to  grow  rapidly,  the 
concepts  of  Rest  Possible  Arithmetic  (RPA)  and  directed  roundings  come  into 
play  (see  references  [1,4]).  RPA  Is  a set  of  algorithms  for  performing  the 
four  basic  arithmetic  operations  while  retaining  the  maximum  possible  amount 
of  accuracy  on  a given  host  machine  (CDC  Cyber  70  Series).  The  directed 
rounding  algorithm  rounds  the  double  precision  result  of  such  an  operation 
into  a single  precision  final  result.  This  rounding  may  take  one  of  four 
directions:  Inwards  (towards  zero),  outwards  (away  from  zero),  upwards  (to- 

wards zero  if  negative,  away  from  zero  if  positive),  or  downwards  (towards 
zero  If  positive,  away  from  zero  If  negative).  The  direction  of  the  rounding 
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is  specified  by  the  computational  envl ronment  In  each  Interval  routine. 

Our  CDC  implementation  of  this  software  package  has  followed  the  sug- 

i 

goat  Ions  of  Yohe  and  Introduced  only  those  changes  necessary  to  accommodate 
the  CDC  software  and  floating-point  arithmetic.  These  changes  are  described 


in  the  UPDATE  deck  on  our  release  tape  (see  Sect  ton  4).  For  the  HPA  routines, 
we  have  Implemented  two  versions.  The  first  version,  called  RPALIR,  contains 
the  following  routines: 

BPAADD  - Yohe ' s addition  algorithm 
RPASUB  - Yohe's  subtraction  algorithm 
BPAMUL  - Yohe's  multiplication  algorithm 
BPAD1V  - Yolie's  division  algorithm 
BPARND  - Yohe's  directed  rounding  algorithm 
All  five  routines  Implement  the  algorithms  given  In  [1)  with  modifications 
made  only  to  allow  tor  the  one's  complement  Cyber  70  floating-point  format. 

In  addition,  we  have  implemented  a second  vet; Ion,  called  RPFt.lR  (resi- 
dent in  RPAl.IB) , which  contains  the  following  routines: 

RPFADD  - Sterbenz'  algorithm  for  an  ssd  addition 

KPFSUB  - Sterbenz'  algorithm  for  an  ssd  subtraction 

BPFMUl.  - Sterbenz'  algorithm  for  an  ssd  multiplication 

RPFDIV  - Sterbenz'  algorithm  for  an  ssd  division 

ID'KRN!)  - a double  to  single  rounding,  similar  to  RPARND 
These  routines  implement  the  algorithms  suggested  by  Sterbenz  [6].  All  these 
routines  are  coded  in  CDC  6000  assembly  language.  In  contrast  to  the  Rl’A 
routines  which  implement  Yohe's  algoilthms  to  perform  best  possible  floating- 
point arithmetic,  the  RPF  routines  use  the  CDC  6000  floating-point  hardware 
to  perform  faster,  reasonably  accurate  computations  that  are  within  one  ulp 
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of  the  corresponding  HI’ A insults;  Initial  numerical  tents  ot  our  Imple- 
mentation Indicate  that  Bl’K  routines  tend  to  require  from  50  to  80  percent 
of  the  execution  time  of  their  Bl’A  counterparts.  The  CDC  loader  may  be 
used  to  subntltute  any  Bl’K  routine  for  anv  Bl’A  routine  If  computation 
time  Is  more  critical  than  computational  accuracy  (tighter  bounds);  the 
LDSET  card  with  SUBS'V  parameter  allows  the  user  to  switch  from  any  Bl’A  to 
anv  Bl’K  routine  (see  Section  h).  I'tu  documented  listings  ot  Hl’Allh  and 
Hl’KI.IB  are  available  from  the  Automatic  Data  Processing  (ADH)  Center  ot  the 
l'.  S.  Armv  Engineer  Waterways  Experiment  Station  (WKS),  1’.  0.  Box  nil. 


Vicksburg,  Miss.  3lM80. 


4.  limitations  of  Interval  Analysis  Package 

The  potential  user  of  this  or  any  other  Interval  analysis  package  should 
take  Into  consideration  the  limitations  of  this  technique.  Some  of  the  limita- 
tions are  Inherent  in  all  Interval  analysis  implementations  and  others  are 
dependent  upon  our  specific  package.  First  of  all,  the  user  should  be  aware 
that  software  interval  analysis  tern's  to  be  very  slow;  cases  run  on  our  CDC 
machine  have  required  as  much  as  100  fold  increase  in  execution  time  over  a 
non-interval  version.  From  10  to  300  fold  increases  of  execution  time  over 
non-interval  versions  have  been  observed  for  various  Implementations  of  the 
Yohe  package  (see  Yohe  (1977)3.  Thus  large  production  codes  involving  signif- 
icant amounts  of  floating-point  computation  are  not  viable  candidates  for  a 
complete  interval  Implementation;  however,  small  portions  of  the  computation 
might  benefit  from  the  use  of  the  package  and/or  the  utilization  of  a hard- 
ware-assisted interval  analysis  such  as  that  currently  being  proposed  by 
Intel  Corporation  for  microcomputers  (interval  versions  would  be  only  2 to  4 
times  slower  than  non-interval  versions'  . 

Another  important  aspect  of  interval  analysis  is  that  the  computed  bounds 
can  be  overly  pessimistic  as  well  as  being  time-consuming  to  determine.  The 
innate  structure  of  Interval  analysis  can  make  its  application  insensitive  to 
situations  in  which  common  mathematical  expressions  (such  as  in  a numerator  and 
denominator)  are  assigned  different  bounds;  this  situation  is  usually  referred 
to  as  the  simultaneity  problem  (see  Sterbenz  (1974)).  The  bounds  for  a par- 
ticular computation  are  dependent  on  both  the  form  in  which  the  computation  is 
expressed  and  the  data  which  is  imput  to  the  expression;  thus  it  is  difficult, 
if  not  impossible,  for  the  user  to  define  in  all  cases  the  "best"  form  of  the 
expression  as  well  as  the  "best"  data  in  order  to  ensure  smallest  growth  In  the 
interval  size  associated  with  the  resultant  computation. 
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Furthermore,  the  tightness  of  the  hounds  Is  dependent  upon  the  arithmetic 
of  the  host  computer  and  possible  numerical  anomalies  Introduced  by  compile- 
time arithmetic. 

The  I/O  routines  provided  by  the  Interval  analysis  package  do  not  easily 
permit  as  much  flexibility  as  that  offered  by  standard  Fortran  I/O  routines; 
also  the  Interval  I/O  la  very  slow  (e.g.  :10  minutes  to  Input  10,100  Interval 
numbers  on  the  CDC  Cyber  72).  This  lack  of  flexibility  can  be  partially  over- 
come by  using  Fortran  I/O  with  explicit  lower  and  upper  bounds  of  Interval 
variables  In  the  argument  list. 

Further  difficulties  can  arise  when  converting  non-interval  arithmetic 
or  logical  IF  sLu.e. cents  to  interval  versions.  Since  the  space  defined  hv 
interval  arithmetic  has  no  total  ordering,  any  of  several  possible  partial 
orderings  could  be  employed.  The  user  must  be  careful  when  selecting  an 
Interval  relational  operator  based  on  a particular  partial  ordering  that  he  does 
not  violate  the  Intent  of  the  original  non-interval  IF  statement.  This  Issue 
is  discussed  in  more  detail  In  Section  S along  with  package  supported  relational 
operators  based  on  two  possible  partial  orderings. 

Another  factor  affecting  the  performance  of  the  Interval  package  Is  the 
accuracy  of  the  Fortran  library  routines.  In  order  to  provide  the  tightest 
possible  bounds  on  these  routines,  a thorough  analysis  must  be  performed  of  the 
accuracy  ot  each  rout lne  across  tin-  out  ire  range  ot  representable  numbers  tor 
each  argument.  Bounds  (not  necessarilv  optimal)  on  the  errors  associated  with 
i'l'C  Fortran  library  routines  used  on  our  host  machine,  provided  hv  Control 
Data , are  available  from  the  ADP  Center,  WKS.  It  Is  this  author's  opinion  that 
a much  more  thorough  analysis  must  he  conducted  hut  such  a task  Is  beyond  the 
scope  of  the  current  project;  it  could  he  done  with  the  aid  of  a multiple  pre- 
cision package  and  with  detailed  knowledge  ot  the  algorithms  and  Implementation 
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features  employed  in  the  Fortran  library  routines. 

Additional  constraints  on  the  use  of  the  Interval  analysis  package  Include 
memory  requirements  and  provisions  for  dealing  with  Intervals  which  are  infinite 
and/or  include  zero.  The  latter  problem  has  been  somewhat  resolved  by  Kahan 
and  others  who  have  advocated  entenslon  of  classical  interval  analysis.  The 
former  problem  caused  some  problems  at  SMU  where  the  maximum  user  memory 
space  is  200gk  words;  even  using  overlays,  the  storage  requirements  of  the 
package  prevented  application  to  medium  or  large  scale  matrix  problems. 


5.  Using  the  Interval  Ana  lysis  I’ackagi 

This  section  focuses  attention  on  the  mechanics  of  running  the  Interval 
analysis  package  on  CDC  computers.  An  overview  of  the  actions  taken  by  the 
control  cards  used  In  this  section  Is  given  In  Section  2.  It  la  assumed  that 
the  potential  user  of  the  package  can  modify  his  Fortran  program  for  use  with 
the  Interval  anal vs  is  package  In  accord  with  the  limitations  discussed  In 
Section  A.  Here  we  will  discuss  the  general  and  specific  form  of  the  job 
control  cards  used,  the  Input  and  output  of  Interval  quantities,  and  an  ex- 


ample run  with  the  given  software. 


The  following  control  card  sequence  Is  typical  for  use  of  'he  AUCMENT 
pi ecomp 1 1 er  : 


* 


* 

a 

* 

* 


* 

* 

* 


JOB;  ACCENT. 

ATTACH,  AUGMENT. 

SKTTL(tim»  limit) 

RFL,  (size  ol  AUGMENT) 

AUGMENT  (,  SF,  SCI) 


(Inst a 1 lat Ion  dependent ' 

(not  required  if  AUGMENT  Is  system-resident; 
FN,  UN,  and  NA  parameters  m.iv  he  required) 

(varies  from  20k  to  200k 
depending  upon  overlay  tactics) 


RFWIND,  SF. 

RETURN,  SCI. 

RETURN,  AUGMENT 
RFL,  50000. 

ETN,  I-SF(,  your  options). 


RFL,  20000. 

ATTACH,  INTUB,  BFAL1B . 

LIBRARY,  1NTLIB,  Bl’AUB. 

LDSKT  (USEP  - BPABKT T/ INTBKPT  (,  your  options)) 
RFL,  (program  sire) 

LC0  (,  FILE  parameters). 


RFL,  20000. 

LIBRARY. 

RETURN,  INTL1B.  BFAUB,  LC0. 


(control  cards  for  post-interval  analysis  processing) 


(end  of  record) 


All  cards  marked  with  * are  optional,  hut  generally  have  the  effect  of  improving 
system  throughput.  Note  that  one  of  the  options  on  the  LOSET  card  could  he  the 
"substitute"  option  to  replace  references  to  BPA  routines  with  references  to 
the  BFF  routines  (see  discussion  In  Sect  Ion  ')  it  the  faster  arithmetic  provided 
by  these  routines  Is  desired. 


IP 


The  entire  job  deck  structure  Is  as  follows: 


(control  card  record) 

(end  of  record) 
description  deck 
♦BEGIN 

program  In  extended  FORTRAN 
♦END 

(end  of  record) 

Input,  If  any,  for  extended  FORTRAN  program 
(end  of  record) 

input,  If  any,  for  post-interval  program 
(end  of  Information) 

At  SMU,  the  software  Interval  analysis  package  is  run  on  a CDC  Cyber  7.’  and 
resides  on  a private  disk  puck  (PRO  10 .16'.  The  specific  deck  setup  for  a run 
with  this  package  Is  as  follows  In  our  environment: 

Bin  number.  name 

ACCOUNT,  ’ character  alphanumeric  account  no,  password. 

SETTL,  appropriate  t im*  limit  in  octal  seconds. 

ATTACH,  AUCMENT/UN-EAUR005,  l’N-UR0101b , M-V.NA. 

ATTACH,  ISTLlB/PN-PROlOlb,  VN-KACROO ' , M-R.NA. 

RFI,,  appropriate  number  ot  words  in  octal. 

AUGMENT  ( INPUT , COMPILE,  CUT  PUT , n -990009,  Al’OSRC) 

RT,  Al'CSCR. 

ATTACH,  REAL  IB/ VN-KACROO 7,  PN-PR01016 . M-R. 

FTN,  T,  PL-9990999. 

LIBRARY,  INTLIB , BI’ALIB . 
l.DSKT  (USFP-RPAbKDT/ 1NTBKDT) 

LC0 . 


♦DESCRIBE  EXTENDED 
AUGMENT  description  deck 
♦BEGIN 

PROGRAM  Name (INPUT,  OUTPUT,  TAPE 5“ INPUT,  TAPEb-OVTPVT , P^RUC-OUTPUT' 


Fortran  source  program 
* END 

h 

9 

data,  if  any 


1 ' 


In  the  above  control  all  command - are  written  beginning  in 


column  1.  the  card  7 represents  a multiple  punch  in  column  l signifying 

% 

end  of  record  and  the  card  b,  represent*  a nil  tiplc  punch  in  column  l in- 

dicatlng  end  of  information.  Note  that  for  the  CDC  PKOC.RAM  card,  a slash  uust  he 
inserted  In  column  l in  order  to  indicate  to  the  AIY.MVN 1 preprocessor  that 
it  is  to  ignore  this  statement  (otherwise  a ,lr  agnostic  error  message  is 
generated!.  The  word  FROGKAM  begins  in  column  7 on  this  card.  If  the  user 
wishes  to  use  the  taster  KPK  intetval  arithmetic  rather  than  the  BPA  interval 
arithmetic  specified  hv  the  above  contvo’  sequence,  then  t wo  changes  ,u  e 
necessary  ' 

replace  ATTACH,  BPAUP.T  N-FACKOO  : , VN-PROlOlb,  N-X. 

with  ATTACH,  BPAM  mi  l A : VN-l  VROO  . VN-l'ROlClb,  NA  . 

and  after  LDSKT (I'SFF-KFAKkV'.  1 \TKklv. ' 

insert  l PSlVySl HS'-KVAA'"  M’VM'i'  '■ AS'.".’.  KWSl'K  KVAMVl  hl'1'’.Vl  Kl'AVlV  Hl'Vl'lV 


The  above  insert  can  be  varied  to  perr  • the  ptogtammet  to  m>.  x the  use  ot  HPA 
and  BFF  arithmetic;  for  example,  the  use  o'  : h.e  control  card 

L11SFT  y SVBST"Bl'AMVL  KlTV.ll  KVAl'lV  - KIT1VV ' would  result  in  KVA  arithmetic 
being  performed  for  addition  and  subtraction  and  KV*'  aiithmetl,  being  pet  termed 
tor  multiplication  and  division. 

From  the  information  provided  hv  the  release  tape  (described  in  Sect  ton  (>' , 
the  programmer  could  set  up  the  package  to  be  u ed  with  the  above  control  se- 
quence employed  at  HVV  ot  could  devise  ant  othet  .ntangement  evte  convenient  to 
the  uset  and  or  his  computet  installation. 

An  Interval  variable  ot  an.iv  is  declare,:  by  the  1NTFKVA1  declaration  which 
when  processed  hv  Al’CMFST  generates  a RKAl  dec lat at  ion  with  an  additional  dl 
mens  loti  introduced  to  provide  space  for  both  endpoints  ot  each  interval  element; 
for  example,  a declaration  of 


INTERVAL  A,  B(3,3) 


would  generate  a declaration  of 

« REAL  A(2) , B(2 , 3 , 3)  , 

where  A(l)  and  A(2)  would  represent  the  left  and  right  endpoints,  respectively, 
of  interval  variable  A,  B(l,l,l)  and  B(2,l,l)  would  represent  the  left  and 
right  endpoints  of  interval  variable  B(l,l),  etc.  Interval  quantities  can  be 
expressed  in  a variety  ot  ways  both  within  the  program  or  a data  input  to  the 
program;  a description  of  these  representations  is  provided  by  Yohe  (1977)  and 
is  reproduced  as  Appendix  C to  this  report. 

Interval  quantities  can  be  input  or  output  via  standard  READ,  WRITE,  or 
PUNCH  statements  and/or  by  the  INTRDF,  INTRD,  or  1NTWR  statements  provided  with 
the  package.  Wien  using  the  standard  Fortran  1/0  routines  with  interval  argu- 
ments, it  is  necessary  to  place  a / in  column  1 (to  tell  AUGMENT  not  to  pro- 
cess this  card)  and  to  refer  to  the  interval  quantities  in  their  expanded  form 
as  generated  by  the  INTERVAL  declaration.  Thus  if  for  the  Interval  quantities 
referred  to  in  the  preceding  paragraph  above  we  wanted  to  output  the  left  and 
right  endpoints  of  A,  B(2,l),  and  R(2,2)  ve  could  use  the  following  statement: 

/ WRITE (6 , 100)  A(1 ,2) , A (2 , 2) , B(l,2,l),  B(2,2,l),  B(i,2,2),  B(!,2,2) 

100  FORMAT  (IN,  6E21-14) 

We  can  proceed  in  a similar  fashion  when  using  READ  or  PUNCH  statements,  i.e. 
slash  in  column  1 and  Interval  variables  expressed  in  their  expanded  dimensional 
form  when  referenced  in  the  argument  list. 

The  format  for  the  I/O  provided  by  the  interval  analysis  package  is  as 
foil ows : 

CALL  INTRDF  (UNIT,  X) 

CALL  INTRD  (UNIT,  FMT,  A,  N) 

CALL  INTWR  (UNIT,  FMT,  A,  N) 


where 


UNIT  • 5(standard  Input)  or  Ofstandard  reload) 

hist  an  laid  print  unit)  or  1 (standard  punch  unit) 

KMT  “ a user  supplied  format  which  must  he  declared  by  the  user  as 

a )-element  Integer  array  for  use  with  Input  statements  and 

as  a 4-element  Integer  array  for  output  statements. 

FMT(l)  • number  of  Interval  variables  In  each  record 

FMTt?)  ■ number  of  characters  to  he  slipped  between 
Interval  variables  ir.  each  record 

FMT(3'  “ numb ei  of  decimal  digits  In  each  interval  variable 

" (for  output^  .'A0+l  b where  Q Is  number  of  decimal 

digits  to  be  displayed  fot  each  Interval  endpoint 

KMT(*'  - (tor  output  only)  Hollerith  printer  carriage  control 

character  ’ ' or  'O' 

For  input,  the  1 NT RDF  routine  obtains  the  next  data  field  from  the  input  stream, 
converts  it  and  stores  the  result  fn  the  declared  Interval  variable  specified  as 
the  second  argument  In  the  call  to  the  routine  (this  variable  is  X in  the  above 
routine  call);  only  one  value  Is  read  by  each  call  to  routine  1NTRDF.  In 
contrast,  INTRP  can  ho  used  to  Input  the  first  N elements  of  an  arrav  of  Interval 
variables  whose  name  is  specified  as  the  third,  argument  in  the  call  to  the 
routine  (A  In  the  above  example';  the  Interval  data  elements  .an  be  expressed 
In  anv  of  the  forms  given  In  Appendix  0.  c numbot  of  elements  to  be  read 
can  be  expressed  explicitly  or  as  any  user  - def tned  integer  variable  given  ns 
the  fourth  argument  in  the  call  to  routine  IMTRD.  For  output,  the  Interpre- 
tation of  the  third  and  fourth  arguments  for  INTVR  Is  essentially  the  same  as 
that  for  I N TRD  except  that  the  Interval  quantities  are  expressed  in  one  form 
in  which  each  endpoint  ts  represented  bv  a l)  decimal  digit  signed  mantissa 
wltb  a three  digit  signed  exponent  as  specified  bv  KMT (3)  above.  Since  our 
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CDC  machine  has  roughly  the  equivalent  of  14  or  15  decimal  digits  In  single 
precision  floating-point  arithmetic,  the  selected  value  for  Q should  probably 
by  less  than  or  equal  to  15.  The  value  of  FMT(3)  should  not  exceed  135  since 
the  width  of  our  printer  line  is  136  (which  Includes  1 character  for  carriage 
control).  Above  we  have  given  a very  brief  account  of  I/O  with  the  Interval 
analysis  package;  extensive  additional  details  are  given  by  Yohe  (1977). 

In  the  remainder  of  this  section  an  example  is  presented  which  illustrates 
the  use  of  the  interval  analysis  package  on  a CPC  Cyber  72  computer  and  how  the 
non-interval  version  of  a Fortran  program  can  be  modified  to  produce  an  inter- 
val version.  The  main  program  SPI.1ST  calls  subroutines  SrLlNT  and  SPLINT 
which  are  used  for  spline  interpolation  problems.  Computational  details  and 
analysis  of  results  using  these  routines  can  be  found  in  the  paper  entitled, 
"Some  Numerical  Computational  Kxperi ence  with  a Software  Interval  Analysis 
Package  and  FORTRAN  Preprocessor  on  Control  Data  Cyber  70  Series  Computers" 
by  Myron  Ginsberg  and  available  as  report  CS  7808  (July  1978)  from  the 
Department  of  Computer  Science  at  Southern  Methodist  University  in  Dallas, 
Texas.  The  non-interval  version  of  the  main  program  and  the  accompanying  sub- 
routines along  with  the  dnyfile  (includes  Job  control  statements  as  well  as 
statistics  about  the  execution  or  the  program’'  are  available  at  UTS,  as  is  art 
interval  version  with  its  dayfile. 

Thi  interested  reader  should  obtain  these  from  WF.S  and  compare  the  non- 
interval and  Interval  versions  to  observe  the  similarities  and  differences. 

It  is  first  noted  that  the  dayfiles  are  substantially  different;  this  occurs 
because  the  non-interval  version  merely  requires  a Fortran  compilation  of  the 
source  program  whereas  the  interval  version  requires,  in  addition,  the  use  of 
the  Fortran  preprocessor  AUGMENT  and  a large  collection  of  interval  arithmetic 
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support  routines.  The  job  control  statements  lor  the  interval  version  are 
in  accord  with  the  discussion  at  the  beginning  of  the  section.  Also  the 
PROGRAM  card  in  the  non-interval  version  has  i blank  in  column  1 whereas  a 
slash  appears  in  column  1 of  the  PROGRAM  card  for  the  interva  version;  the 
slash  indicates  to  AUGMENT  that  it  is  not  to  attempt  to  alter  this  card  but 
instead  to  pass  it  on  intact  (but  without  the  slash)  to  the  Fortran  compiler. 
Similarly,  other  statements  in  the  interval  version  with  a slash  ir.  column  3 
(for  example,  the  PUNCH  statements)  are  not  to  be  processed  by  AUGMENT.  In 
the  earlier  part  of  this  section  c o’  slashes  was  mentioned  wifh  I/O  state- 
ments such  as  REAP  or  WRITE  when  interval  variables  are  present  in  the  argument 


Observe  in  the  interval  version  ’hat  there  are  INTERVAL  statements 

present  at  the  beginning  of  the  program  and  each  subroutine  ar.d  that  all 
single  percision  f loat lag-point  variables  used  in  each  routine  are  declared  ir. 
the  accompany ir.g  INTERVAL  statement.  Such  statements  when  processed  by 
AUGMENT  generate  a REAL,  statement  in  which  each  variable  is  changed  into  a 
variable  of  the  same  name  but  with  an  additional  dimension  added  as  the  new 
first  dimension  of  the  original  variable;  this  action  is  necessary  to  provide 
storage  for  the  lower  and  upper  bounds  associated  with  each  interval  quantity. 
.Tie  REAL  statements  generated  by  AUGMENT  for  our  example  under  discussion  are 
also  available  from  the  AIM’  Center,  WHS. 

Another  alteration  of  a non-interval  program  required  to  transform  it 
into  an  interval  version  involves  the  re-writing  of  arithmetic  and  logical  IF 
statements  which  employ  interval  quantities.  Since  interval  arithmetic  has  no 
total  ordering,  any  of  several  possible  partial  orderings  could  be  used  to  pro- 
vide definitions  of  relational  operators.  Two  sets  of  such  operators  are  pro- 


vided  by  the  interval  analysis  package  - value  operators  and  sat  operators. 
The  value  operator  definitions  are  as  follows: 


.VE. 

where 

[a  ,b  ] 

as 

[C.d] 

if 

and 

only 

if 

a = b “ c * d 

.VNE. 

where 

[a,b] 

[c.d] 

if 

and 

only 

if 

they  are  disjoint 

• VLT. 

where 

[a,b] 

< 

[c,d] 

if 

and 

only 

if 

b < 

c 

.VLE. 

where 

[a.b] 

< 

[c.d] 

if 

and 

only 

if 

b < 

c 

.VGT. 

where 

[a.b] 

> 

[c.d] 

if 

and 

only 

if 

a > 

d 

. VGE. 

where 

[a,b] 

> 

[c.d] 

if 

and 

only 

if 

a > 

d 

The  set 

operator 

definitions  are  as  follows: 

.SEQ. 

where 

[a.b] 

= 

[c.d] 

if 

and 

only 

if 

a = 

c and  b = d 

.SNE. 

where 

[a,b] 

4 

[c.d] 

if 

and 

only 

if 

a 4 

c or  b 4 d 

. SLT. 

where 

[a,b] 

< 

[c.d] 

if 

and 

only 

if 

a < 

c 

.SLE. 

where 

[a,b] 

< 

[c.d] 

if 

and 

only 

if 

a < 

c 

.SGT. 

where 

[a,b] 

> 

[c.d] 

if 

and 

only 

if 

b > 

d 

. SGE. 

where 

[a,b] 

> 

[c.d] 

if 

and 

only 

if 

b > 

d 

The  reader  should  note  that  for  the  value  operators , if  two  intervals  are  non- 
empty and  with  corresponding  endpoints  equal,  then  neither  equality  nor  non- 
equality of  these  intervals  is  true;  in  such  cases  the  user  could  decide  to  use 
the  set  relational  operations. 

Several  alterations  of  the  non-interval  versions  of  IF  statements  are 
present  in  the  program  under  discussion  (e.g.  see  lines  9510  or  9590  in  routine 
SPLINE).  For  example  line  1580  of  routine  SPLINT  in  the  non-interval  version 
is  transformed  from 

IF (XP  - X(l))  100,  170,  110 
to 

IF (XP  .VGT.  X(l) )G0  T0  110 
IF (XP  .SEQ.  X(1))G0  T0  170 

in  the  interval  version.  The  user  must  alter  all  logical  and  arithmetic  IF 
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statement  s Involving  inleivul  variables.  When  the  programmer  Is  making 
.such  changes  lie  must  be  e.irelul  t o ric*  i « . I tin*  .'j>i  late  operatoi  which 

best  let  leets  tin*  intent  et  t be  eeoi|  at  *n  in  tin  inut  - ini  er  val  version  of 
tin*  IK  statement  in  the  eentext  ot  t b<*  program;  stieli  soled  Ion  ot  Interval 
operators  Is  particularly  crucial  tot  IF  tt.i  foments  Involving  terminal  Ion 
anil /or  tolerance  criteria  when  an  luappt  opt  i at  e operator  selection  conlil 
stgnl f leant Iv  alter  the  behavior  of  tbe  Interval  version  of  'ite  program  from 
its  non-interval  i ounterpar t . 

As  discussed  earlier  In  this  section,  tbe  user  has  several  ways  to  pro- 


vi ilo  I/O 

f or 

i nc  e 

rva  1 

quant 

It  l ’ 

and  a i cuph* 

ol  t ln*fu*  art*  pi  fSt’iu  In  our 

example  p 

*'<>»**■ 

am . 

Note 

• i hat 

an  v 

>i  d .ii.ii  v PK  i ST 

st  at  omen : is  ust ‘i!  i t ht*  ma  Iti 

prop, rim  Ik'  output  the  non-  interval  quality  VI’  whereas  in  the  interval  version, 
r out  Inc  1 N rVJR  item  tin*  Interval  ana  1 v s l ■ package  l ■ . used  to  on  t pu  t the  1 n - 
terval  variable  \T.  Also  not  o ’hat  the  format  iucJ  with  rout  Ini'  l N'TWR  is 
provtilod  hv  integer  arrav  KMT  whieli  is  explieltlv  deelarei!  and  Initialised 
hv  the  user  in  t he  rout  lire  in  which  INTu'K  i •;  Invoked.  \lso  observe  that 
Interval  varlahW’  VP  is  output  hv  the  use  ot  a I'l'Nr'H  *;t  c ’ ement  hut  that  in 
tills  statement  a slash  appears  in  column  ! Al'eMl'T  will  not  pieces*  the 

statement'  and  that  the  lower  and  upper  hounds  of  YV  (name',  v VP(1>  and  YP(2)  , 
respectively)  are  explicitly  referenced  in  the  argument  list. 

In  completing  the  transformation  of  the  non- Intet val  program  Into  an  tn- 
terval  program,  the  usei  can  utilize  several  operators  that  are  supported  hy  the 
Interval  analysis  package  and  ippticah'e  to  Inteiva!  quant  1 1 t os . A list  ot  all 
such  opei  a t or  s and  t lit*  1 1 del  luit  tons  i s avaiable  at  Wf  S . In  our  sample  program 
we  have  used  MIX’  V and  I l.  01  t o provide  the  midpoint  and  length,  rospeetlve.lv,  ot 
any  Interval  quantity. 
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From  our  discussion  above  the  user  should  lealize  that  to  change  his  non- 
1 interval  program  Into  an  interval  version  on  CPC  machines,  he  muse  perform  the 
following  steps: 

1.  Use  the  Job  control  statements  (or  equivalent)  to  invoke  the 
AUGMENT  preprocessor,  interval  analysis  package,  and  Fortran 
coropt ler . 

2.  Place  a slash  in  column  l of  the  PROGRAM  statement. 

3.  Introduce  an  INTERVAL  declaration  in  each  routine  for 

all  single  precision  floating-point  variables  (and  arrays) 
appealing  in  that  routine. 

A.  Change  all  arithmetic  and  logical  IF  statements  Involving 
interval  quantities,  replacing  Fortran  relational  operators 
with  the  appropriate  interva'  value  or  set  operators. 

3.  For  1/0  interval  quantities  use  the  interval  supported 

routines  (1NTRPK,  INTRP,  INTU'K)  ot  the  regular  Fortran  I/O 
routine  (hut  with  slash  in  column  1 and  explicit  reference 
to  lower  and/or  upper  bounds  oi  each  interval  quantitv  In  the 
argument  list). 

6.  Introduce  any  desired  interval  oriented  operators  supported  bv 
the  interval  package. 

7.  Select  an  appropriate  time  limit  (use  SETT1,  instruction)  and 
field  length  (use  RF1,  Instruction1'  tor  the  interval  version  of 
your  program,  keeping  in  mind  that  the  interval  version  may  re- 
quire from  50  to  100  times  (or  more  the  execution  time  of  the 
non-interval  version  and  that  the  storage  requirement  for  all 
Interval  variables  and  arrays  Is  doubled  from  that  ot  t lie 
non-interval  represent  at  ions. 

After  the  non-interval  program  has  been  converted  into  an  interval  source 
program  using  steps  1 through  6 above,  the  AUGMENT  preprocessor  generates  another 
Fortran  program  which  invokes  calls  to  subroutines  provided  by  the  Interval 
analysis  package  in  order  to  execute  the  interval  program. 


b . Format  of  the  Release  Tape  for  Distribution  of  the  CPC  Version 


The  interval  analysis  software  is  currently  contained  on  two  CDC  UPDATE 
libraries.  One  library  contains  the  BPA  primitives  and  the  other  library 
contains  all  other  routines  Including  the  AUGMENT  routines.  The  format  of 
the  release  tape  is  as  follows: 

7-track 
556  BP  1 


K.RONOS  Internal  Format  (512  word  blocks) 
unlabeled 


File  No . 


Dear r ipt ion 


1 BI’A  and  Bl'F  Primatives  (Source  Card  Images) 

2 INTERVAL  and  AUGMENT  (Source  Card  Images) 

3 BPA  and  BPF  Primatives  (UPDATE  Source  Format) 

4 INTERVAL  and  AUGMENT  Package  (UPDATE  Source  Format) 
Files  1 and  2 are  card  images  which  can  be  read  by  almost  any  machine.  Files 

3 and  4 are  to  be  used  with  UPDATE  to  create  new  UPDATE  libraries.  They  are 
basically  the  same  as  files  1 and  2 except  they  have  UPDATE  directives  embedded 


within  them. 


FVAl.UATtON  OF  INTFKVAI  ARITHM1 TH' 
FOR  THK  HU'  CYUKR  70  SYSTEM 


Problem : To  solve  for  x the  linear  algebraic  system  of  equations 

Ax  " b 

where  A is  full  (turn-sparse)  matrix.  Test  data  for  A Is 
100  x 100  ease  from  WFS 


Code  used  and  references: 

FORTRAN  code  from  WES  called  GAUSS  and  wt  It  ten  hv 
F.  T.  Tracy . 


Conditions  under  which  test  was  conducted: 

The  F'l’N  A. 6 FORTRAN  compiler  was  used  with  no  optimi 
ration  (i.e.  opt  * 0).  Noninterval  and  interval  (both 
HP A and  RPF)  cases  were  run.  T imtnR  figures  were  ob- 
tained by  using  the  SFi  ON!'  function  to  measure  the 
CPU  time  for  each  run  excluding  time  (or  I/O.  The  INTRO 
and  1NTWR  routines  were  toted  with  all  variables  declaied 
Interval.  Inequalities  involving  interval  variables 
used  the  interval  operators  VI. I and  VCF.  to  replace  the 
noninterval  operators  IF  and  OF,  respectively. 


Execution  time  for  noninterval  version: 


14.763  seconds  CPU  time  (excludes  time  involved 
with  I/O) 

Execution  time  tor  BPF  interval  version: 

1139.251  seconds  CrU  time  (18.99  minutes]  (ex- 
cludes time  Involved  with  I/O) 


Execution  time  lor  BPA  interval  version: 

1424.297  seconds  CPU  time  (23.74  minutes]  (ex- 
cludes time  involved  with  I/O) 


Observat ions : 


1.  The  BPF  interval  version  requires  80%  of  the  execution  time 
used  by  the  BPA  version  and  produces  competitive  bounds. 
(This  compares  with  81%  for  DECjlMP  and  SOl-VE  code  with  WES 
data . ) 

2.  The  HPA  and  BPF  Interval  versions  take  96.477  and  77.169 
as  much  CPU  time,  respectively  as  the  noninterval  version. 
(This  compares  with  83.863  and  68.072  for  DEC0MP  and  S(/1,VK 
code  on  this  same  WES  data.) 


3.  The  BPF  Intervals  for  the  x elements  are  wider  than  the 
corresponding  BPA  intervals;  79%  of  these  intervals  have 
a difference  in  width  with  corresponding  BPA  intervals  on 
the  order  of  10_t|  to  10"^*.  This  is  the  same  distribution 
of  interval  width  differences  as  occurred  with  the  same 
problem  with  the  DECO MV  and  Sl/LVE  code,  although  the  specific 
values  of  the  Interval  width  differences  vary  slightly 
from  those  of  the  DECAMP  and  SOhVK  code.  The  distribution 
of  BPF  - BPA  interval  widths  for  x elements  is  as  follows: 


difference  in  width 
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Observations:  (cont.) 
4 


The  distribution  of  BPA  interval  widths  for  the  components 
of  the  solution  vector  x is  the  same  as  for  the  BPA  intervals 
using  the  same  data  and  the  DEC0MP  AND  S0LVE  code.  The  dis- 
tribution is  as  follows: 


difference  in  width 


60% 
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The  distribution  of  BPF  interval  widths  for  the  components 
of  the  solution  vector  x Is  the  same  as  for  the  BPA  interval 
widths  as  well  as  for  the  BPA  and  BPF  interval  widths  from 


60% 
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same  data 

used  with  the  DECAMP  and  Sifft.VE  code: 
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Since  the  distribution  of  widths  of  intervals  is  the  same  for 
BPA  and  BPF  versions,  it  is  therefore  justified  to  use  BPF 
bounds  especially  since  BPF  bounds  require  less  execution  time 
than  BPA  bounds  (80%  of  BPA  execution  time  for  this  problem). 


The  computational  results  (x)  (both  interval  and  noninterval 
versions)  for  both  the  GAUSS  code  and  the  DF.C0MP  and  S0LVE 
code  are  remarkably  similar  which  is  somewhat  reassuring  con- 
sidering two  different  codes  (with  slightly  different  algorithms) 
were  used  on  a problem  that  seems  to  slightly  ill-conditioned 
(determinant  \ .361  x 10~20  and  condition  number  ^ 8176).  Pos- 
sibly with  a better  conditioned  A matrix,  the  bounds  with  both 
codes  might  have  been  tighter  than  the  order  of  magnitude 
distribution  indicated  in  observations  4 and  5. 


The  computed  results  for  x , the  residuals,  interval  widths  and  midpoints 
of  x components,  ns  well  as  CPU  timings  are  given  for  the  noninterval 
and  interval  (BPA  and  BPF)  versions  in  the  attached  listings. 
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Execution  time  for  noninterval  version: 


18.902  seconds  CPU  time  (Includes  calls  to  DECAMP 
and  S01.VE) 


Execution  time  tot  BPK  Interval  version: 

1286.701  seconds  [.'1.445  minutes]  CPU  time  (Includes 
calls  to  DEC0MP  aid  S0W) 


Execution  tine  for  BP A interval  version: 

1586.217  seconds  [26.47  minutes]  CPU  time  (includes 
calls  to  PEC0MP  and  S0LVF) 


Obset vat  ions : 

1. 


The  HPF  Interval  version  requires  81?  of  the  exe- 
cution time  used  by  the  BPA  version  and  produces 
compet i t ive  bounds . 


The  BPA  and  BPV  interval  versions  take  83.865  and 
t>8.072  as  much  tine,  respectively,  .is  the  non- 
interval version. 


The  BPE  intervals  for  the  x elements  are  wider  than 
the  corresponding  BPA  intervals;  79*  of  these  in- 
tervals have  a difference  in  width  with  corresponding 
BPA  Intervals  on  the  order  of  10  1 to  !0'H.  The 
distribution  of  BPK  - BVA  Interval  widths  is  as 
fol lows : 
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Observations : (cont.) 


4.  Tlio  distribution  of  BPA  interval  widths  for  the 
components  of  the  solution  vector  x is  as  follows: 
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. The  distribution  of  BPF  interval  widths  for  the 
components  of  the  solution  vector  x is  as  follows: 


iO 

10 

10 


-5 


-6 


10 

10 

10 

10 


-8 


-9 


-10 


X 

XXXXX  XXXXX  XX XXX  X 

XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX 

XXXXX  XXXXX  XXXXX  XXXXX  XXXXX 

XXXXX  XXXXX  XXXXX  X 

XX  XX 

XXX 


We  observe  that  the  distritution  of  widths  of  intervals 
is  the  same  for  BPA  and  BPF  versions,  therefore  further 
Justifying  the  use  of  BPF  bounds  since  the  BFF  bounds 
also  require  less  execution  time  than  BPA  bounds  (81T. 
of  HFA  execution  time  for  this  problem). 


6.  We  observe  that  the  computed  value  of  the  determinant 
of  A is  small  (v.Jbl  x 10--^)  and  the  extlmate  of  the 
condition  number  (x  817b)  Is  fairly  large,  thus  indi- 
cating some  ill-conditioning  in  the  problem. 


The  computed  results  for  x,  the  residuals,  interval  widths  and  mid- 
points of  x components,  as  well  as  CPU  timings  are  given  for  the 
noninterval  and  interval  (BPA  and  BPF!  versions  ft:  the  attached 
1 is  t i ngs . 
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'J.  Problem: 


To  solve  for  x the  linear  algebraic  system  of  equations 


Ax  “ b 


where  A Is  sparse. 


Test  data  for  A Is  10  x 10  block  trldlagonal  matrix  with 
10  x 10  blocks  where 


C 1) 

R C 1U 


-1  4 


cr-p 


-1.5 

0 - -1.5 

O 


1.5  J 


Code  used  and 
references : 


NSl'lV  - a fortran  routine  which  uses  sparse  Gaussian  elimination 
with  partial  pivoting  and  war.  written  bv  Andrew  H.  Sheruiau 
(University  of  Texas-Aust  tn> . The  code  is  listed  in  NST’.V:  A 
FORTRAN  Subroutine  lot  Spat  se  Gun  .s  ( an  K imt-aat  ion  with  1'artlal 
Pivoting.  Report  TR-65 , Department  of  Computet  Sciences,  University 
of  lexas . Austin,  Texas,  January  ! l>  ' / . Also  see  Algorithms  tot 
Sparse  Gaussian  rumination  with  Partial  Pivoting,  Report  No. 
UIUCUCS-R- 76-81 7 , Department  of  Computer  Science,  University  oi 
Illinois,  i'rbana , Illinois,  I1*?!'.  both  reports  were  written  bv 
Andrew  II.  Shot  man. 


Conditions  under  which  The  FTN  4.b  FORTRAN  compiler  was  used  with  no  optimisation 

test  was  conducted:  (l.e.  opt  * 0).  Nonluterval  and  Interval  (both  BPA  and 

Bl’F)  cases  were  run.  Timing  t igures  were  obtained  by  u'  lg 
the  SECi)WD  function  to  measure  the  CPV  time  for  each  ru. 
excluding  time  for  1/0.  The  1NTWK  routine  was  used  with 
variables  declared  Interval  In  IV  statements  Involving 
interval  quantities,  the  operators  Et),  CT,  and  LT  are  re- 
placed by  SEQ,  VOT,  and  V1.T , respectively. 

Execution  time  for  .34  seconds  CPU  time 

noninterv.rl  vers  ion: 

Execution  time  for  V 1 . hbT  seconds  CPU  time 

RPF inter  va  ’.  vor s Ion  : 

Execution  tine  fo_r  1 3 . 4 p 3 seconds  CPU  time 

BPA  tut  e vv a ver  s ion  : 

Observations:  1.  Slow  interval  version  (BPA'  requires  41.0?  times 

as  much  CPI  time  as  noninterva'  version. 

2.  Fast  interval  version  (BPF)  requires  34.31  times 
as  much  CPV  time  as  noninterva!  version. 

3.  The  fast  interva1  version  (BPVi  requires  B3.S4?  or 
slow  (.BPA'  Interval  Cl'V  time. 

4.  The  ratios  for  fast  and  slow  interval  versions  to 
the  noninterval  version  are  smaller  tot-  this  100  x 100 
sparse  case  than  lor  100  x 100  tronsparso  VKS  case  us'  g 
DEC(?MP  and  SOLVE  (in  that  situation  BPA  version  war 
S-.  times  as  slew  as  nonintot val  version  and  BPF  version 
is  oP  times  as  slow  as  noninterval  version).  As  a 
matter  of  tact,  comparing  tire  specific  nonsparse  and 
sparse  100  x 100  cases  mentioned  here,  we  note  that 
sparse  ratios  are  1/2  those  of  the  cort esponding  non- 
sparse  ratios. 

H 

< 


10.  Prob  \ om : 

To  applv  interval  analysis  to  the  solution  of  a linear  system  of 
equations.  Ax  - b,  wheio  A is  non  singular  but  nearly  singular. 

Wo  want  to  observe  for  slight  perturbations  of  b,  the  effects  on 
the  sice  and/or  midpoint  of  the  generated  Intervals  associated 
with  the  resulting  solution  and  the  relationship  between  these 
values  and  those  corresponding  to  the  original  problem. 

Code  used  and  references: 

DKO>M.P  and  S01.VV  from  Computer  Methods  for  Mathematical  Compu- 
tations by  George  E . Forsythe , HlcHeel  k.  Malcolm,  d Clave  B. 

Me ler , Prentice-Hall , Fng’.ewood  Cliffs,  New  Jevsev , 1911,  pp. 

52-55. 

Cond  '■  i or.s  under  which  tes:  was  conducted  : 

The  VTN  4.0  FORTRAN  compiler  was  used  with  no  opt  imlr.it  ion  (i.e. 
opt  * 0).  Non -interval  and  interval  (both  RPA  and  i'lM'N  cases 
were  run.  Timing  figures  were  obtained.  b\  using  the  SV.Cv'N!)  func- 
tion to  measure  the  CP’d  tire  for  each  run  excluding  time  for  1/0. 

"‘•o  IN  :'K''  IWTWR  routines  were  used  with  all  variables  declared 
Interval.  (equalities  involving  interval  variables  used  the  in- 
terval operators  SV'Q,  VC",  Vl.T,  and  Vl.F  to  replace  the  non-interval 
operators  FQ,  OP,  IV,  LV , respectively. 

Note:  The  test  problem  used  in  this  experiment  is  defined  hv  O.J.  Tee 

in  an  article  entitled,  "A  Simple  Example  of  an  1 11 -Cond it lone. 
Matrix,"  S1CNVM  Newsletter,  Volume  7,  Number  3,  October  1°?/, 
pp.  19-70. 


The  suggested  test  matrix  A is  of  order  1 and  has  a determinant 
of  1 

11  10  Is  ' 


A - 


11  -13 

13  -t>b 


1 ct 


1.001 

0.999 

1.001  j 


then  the  exact  solution  of  Ax  - b is 
I -0.683  'J 
x - ' 0.843 

L 0.006  J 


bet 


b -t  db 


equation  for  this  experiment  is  Ar  » 


fb 


then  the  perturbed 


b + 3b. 


Comput.it  Iona  I res  ul  ts : 

Noninterval  Version*  for  Ax 
Computed  values  for  x 


opt-0 
truncated 
ar  1 thinet  Ic 


f - . 6830000000364  58 
. 84  30000000399 1 1 
V .00600000000012932 


opt"0 

rounded 

ar ichmetic 


( -.682999999999616 
!j  .84299999999957' 

's  .00589999999999803 


b 


Noninterval  Version*  for  Az 


Computed  values  for  x 


opt=0 

truncated 

arithmetic 


1.00000000005106 
\ -1.00000000005591 

^ -1.S11683976283S3  X 10~l* 


b + 6b 


opt”0 

rounded 

arithmetic 


r ,999999999996479 
( -.999999999996138 

V .133226762955049  X 10 


Computational  value  of  determinant  of  A = . 999999999999545, 
Estimate  ot  condition  number  of  A * .164135602084493  X lO*'. 


*Noninterval  versions  were  also  run  under  opt  « 1,2,  and  2 Uo  with 
rounded  and  truncated  arithmetic;  truncated  results  under  all  option 
gave  exact  agreement  as  did  rounded  results  under  all  options. 


Computational  results:  (cont.) 


. ' 


Interval  Version*  for  \x  = b 

BPA  computed  intervals  for  x 

(- . 68000000274032 ,-.682999999775020) 

(.842999999762525, .843000000291192) 

( . 00599999999922828 , . 00600000000094667) 

BPF  computed  intervals  for  >c 
(-.683000000277098, -.682999999773805) 

(.842999999761271, .84300000029^467) 
(.00599999999922420,-00600000000095733) 

BPA  midpoint  of  x length  of  x 

-.683000000024528  .499010610610640  X 10~9 

.843000000026858  .528665111687587  X 10_ 

.600000000008749  X 10~2  .161837544193290  X 10~ 

BPF  midpoint  of  x length  of  x 

-.683000000025451  .503291630593594  X 10_9 

.843000000027871  .533194821628058  X 10~9 

.600000000009077  X 10-2  .173311365259110  X 1C-11 

Total  CPU  time  for  BPA  version  ».198  seconds 

Total  CPU  time  for  BPF  version  ~.164  seconds 

Total  CPU  time  for  Noninterval  version  .005  seconds 


I ’ 

i.i 

4 

4 


1 

4 


*Interval  versions  were  run  for  BPA  and  BPF  under  opt=0,l,2,  and  2 U 0 
with  rounded  and  truncated  arithmetic.  Truncated  and  rounded  interval 
bounds  were  identical  under  all  options  for  BPA  and  likewise  for  all 
BPF  versions. 


Computational  results:  (coot.) 


Interval  Version*  for  A z * b + fib 
BPA  computed  intervals  for  z 
( . 999999999947554 , 1 . 0000000001 I9t>9) 

(- 1.000000000 12800, 99999999994 5636) 

(-••41  '>667  500471075  X 10~ 1 2 1 75859327122378  X 10~12) 

BPF  computed  intervals  for  7. 

( . 998999999944471  > 1.00000000012126) 

(-1.0000000001 2963, -.99999999994 2 139) 
(-.420996570989935  X 10"! ‘ 18651 7468160098  X 10~12) 


midpoint  of  ;• 

lent, th  of  7 

1.00000000003362 

.172128977737884  X 10-9 

-1.00000000001681 

.182353637705472  X !0~g 

. 1 1990408h674J49  } 

c uf12 

.59126827593453  X 10'12 

midpoint  of  z 

length  of  7 

1. 00000000003286 

,17678303]657113  x 10~9 

-1.00000000003598 

.187284854291647  X 10-9 

- .117239551414919  > 

i io~12 

.617514039150033  X 10~V2 

Total  CPI'  time  for  BPA  version  .197  seconds 
Total  CPU  time  for  BPF  version  - .163  seconds 
Total  CPI'  time  for  noninterval  version  ~ .005  seconds 


•Interval  versions  were  run  for  BPA  and  BPF  under  opt"0,l,2,  and  ?C,0 
with  rounded  and  truncated  arithmetic;  truncated  and  rounded  interval 
bounds  are  identical  under  all  options  for  BPA  and  likewise  for  all  APF 
vers  Ions . 


18 


. • .. 


Obsoi vat  ions : 


BPF  ( I’ll  time 
HIT  CIT  t Lino 


<•  ft  It 


(This  compares  with  8l»  loi  the  samo  emit*  when  used  with 
100  X 100  WV.S  case.) 


BIT  CPU  t (me  _ , , 

Noninfeival  Cl’U  t tmo 

(Thi:-  compares  with  f S lot  the  same  code  when  u-oj  with 
100  X 100  WKS  case.) 


Of . s 

Non < nt  o v.i  ! i IT  t (mo 

(This  compares  with  I'1.'1  let  (lie  same  code  when  teed  with 
100  X It’d  Vi's  case.-) 


. The  tat  ivs  ot  the  widths  el  the  intervals  . ort  ospondtup,  to 
the  ,\  and  .’  solutions  at.  appi  ox  ima  t e 1 v the  sane  toi  all 
components  for  BOA  and  Hl’V  versions,  t.e. 


width  of  \ . Intel va  1 
w id  t h ot  i nt  i i va  1 


Thus  the  inlet  vats  it.soeiated  with  the  perturbed  problem 
are  aetuallv  sm.tllei  than  the  cot  respond  i up,  intervals  as  so  • 
c t a t ed  i.  . t h the  x se  1 u t I on  and  the  .'fere  do  not  aha  i the 
us ei  . Although  the  l if, lit  hand  side  ot  A.’-lvKMi  has  been  pet 
Imbed  only  .1’’  t t imi  that  ef  h in  Ac  h , lit.'  s.'lul  ion  . is 
aetuallv  pci  Imbed  '.V0"-  (see  Ot>set  vat  ten  S).  Thus,  in- 

terval widths  in  this  expel  iiaeut  do  not  w.itn  the  uset  ot  the 
.11  condiltonin.  ot  thi  ptoblem  luole  that  the  estir.iate  of 
the  coiulii  ion  turn  .a  .-1  A • It's,  lib  which  indicate'  the 
problem  Is  ill  coti.l  1 1 i one.O  . As  a mat  lot  oi  tael,  the  n- 
tetv.il  leapt’  s 'ei  \ and  computed  ceneoneuts  ate  v.a  \ small 
(on  the  oi  del  o1  10  IJ  to  i 0 1 •)  and  tliet  .h\  do  not  alert  the 
liset  of  nun  el  f a l di  I I , 111  v . 


The  .ihso'.iiio  value  ol  the  telafive  e\  rot  ot  the  midpoints  oi 
the  r intervals  .tt  e bet  let  Indicators  of  the  unmet  leal  I lit. 
tti.it  I on , i . e . 


Yxi 


•i  Is. 


t I 
\ , 


I 


i Me 


1 . 00 


oi 


Observations:  (cent.) 


We  note  particularly  for  the  larger  components  of  x,  the 
absolute  value  of  the  relative  error  In  the  midpoint  of  r. 

Is  on  the  order  of  >218%  (i.e.  for  first  z component ^ 246X 

and  £01  second  z component  219*)  and  this  is  caused  by  a 
.1%  variation  in  the  right  hand  side  of  the  original  system 
Ax”b.  Thus  the  midpoints  of  the  z component  intervals  are 
better  indicators  of  the  numerical  difficulty  than  the 
widths  of  these  intervals. 

6.  Examination  of  the  absolute  value  of  the  relative  error  for 
the  z values  in  the  noninterval  version  of  this  problem  in 
dtcatc.s  approximately  the  sraae  values  given  in  observation  B 
above  when  dealing  with  midpoints  of  the  generated  solution 
intervals.  Thus  the  tionioterval  solution  was  .furl  as  good 
an  Indicator  of  numerical  difficulty  as  wore  the  midpoints 
of  the  interval  solution  and  has  the  added  advantage  of  re- 
quiring less  time  and  storage. 

7.  Note  '.hat  the  residuals  associated  with  the  z solution  (see 
attached  listings)  indicate  values  on  the  order  of  10~^ 
which  might  deceive  the  user  Into  hellev'ng  that  the  results 
are  verv  accurate.  The  user  should  realize  that  for  ill- 
conditioned  problems  such  as  this  one  (tecall  that  the  esti- 
mate ol  the  condition  number  is  164,116),  small  residuals 

do  not  indicate  that  the  computed  solution  Is  close  to  the 
exact  solution. 

8.  The  interval  solution  computed  for  the  perturbed  problem 
(l.o.  Az“b+6h)  does  not  warn  the  user  of  numerical  diffi- 
culty mainly  because  the  widths  of  the  resultant  solution 
bounds  are  dependent  upon  the  number,  tvpe.  and  sequence 
of  arithmetic  operations  involved  in  the  linear  equation 
solver  but  do  not  take  into  account  the  condit ion  number  of 
the  amatrix  in  the  problem. 

9.  The  attached  listings  for  the  noninterval,  RPA,  and  BIT 
versions  of  this  test  for  the  original  and  perturbed  problems 
give  the  computed  solution,  timings,  and  residuals. 
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11.  Problem;  To  apply  t ho  WES  SPI.INF  Inli'i  pol.ti  I.mi  routine  (in  both  non  inletval 
nn<l  interval  versions)  to  UTS  4 point  and  '0  point  data. 

lode  used  and  references: 

The  WES  SPLINE  routine  is  adapted  from  work  by  Grevillo  (see 
TS  Report  No.  89),  U.S.  Arm-  Mathematic  Research  Center.  University 
of  Wisconsin,  June  19<>8).  The  documentation  for  the  WES  SPI.TNE 
routine  Is  contained  in  WE''  1 : is.  e 1 laneons  Repot  t No.  0 71  2 
(July  which  is-  available  from  the  Computer  Analysis  Branch, 

AUP  Center,  U.S.  Atn.v  Wateiw.vs  Experiment  St  it  ion,  Vicksburg, 
Mississippi.  The  Wl'>  routine  SPi  1 NY  written  hv  Jav  Cheek  et  W1S 
is  al  ;o  used  in  this  ex  'et  intent  o compute  V and  V’  values  for  riven 
\ values  based  on  the  spl  iru  p.onot  sled  bv  routine  SPUN!  '.be  ‘PI. INI 
routine  is  discussed  ie  tbo  aforementioned  W't  S paper  at  .'  document.! 
tlou  (s  given  in  the  att.n  lied  1 1st  ing. 

Condition-  undot  which  the  test  was  conducted' 

The  KTN  4.(t  EylUT'tAN  compll-'t  was  u-  ed  with  no  opt  Inti. "at  ion  (i.e. 
opt  ••  «>) . Both  Bl'A  and  BE  interval  vet-ion--  i et  - run  to  compare 
with  the  nonlnt  or\  a 1 ca  -t  . Interval  vat  table--  were  output  r • up, 
the  TWTVK  rimt  ittc.  Timim-  t • p.ures  were  obtained  h\  usim  the 
SF‘  ^N'l)  function  to  nte.isui  i the  C”U  time  Irequal  \ ol  v i up. 

Interval  variables  u--cd  the  Interval  operate--  *UT,',  VcT,  .tin  VI. I 
to  replace  the  noninttrval  opt  ttors  and  ' . t c- pi-  tlvclv. 

KES  provided  the  -i  -point  and  10  pc  nt  data  u-ed.  w -.  : h the  ' c t . 

SP1.TNT  was  used  to  produce  A and  V'  va ues  for  cortesp.-m  ' np,  \ 
values  spaced  at  in.  foments  if  . ."j  aero-.-,  tb.  i - t.  - it-  - -.>*«•<*( 

fled  bv  tin'  4-point  and  10-point  data.  P-.e  ••  -point  d.  ta  pan-- 
X"1.6  1 1'  10  and  the  10-polnt  data  span:-  X 0 to  .’t'.  For  X value*: 
oat  s -E-  the  tattle  t'i  tin'  >.lven  data  set,  routine  S IV  VI  extrapolates 
linoar.lv  from  the  I list  et  last  data  point  using  the  .lop.'  of  the 
spline  at  that  endpoint. 


Comout.it  Iona  1 t osu  Its: 

Cri!  timing  (In  seconds)  tor  SPLt'E: 
With  Function  Call  overhead 
Nonintetval  .007 


4-pet  tit  Case 
Without  Function  Call 
.007 


('vet  bead 


10 ’A 


lil'F 


. 4 : l 


.46  7 


. TO) 


CPU  timing  (in  second;:  1 tot  SPI.INF: 
With  Function  Call  Overhead 


1 0-  po  I nt  Cas.  o 
W 1 1 hout  Eunct  i on  Ca 1 1 


Overhead 


Nonintetval  . 070  .014 

Bl'A  2.690  2.684 

BPF  2.161  2.1S7 

Add  i'  I ona  1 comput.'tt  ional  tesults  tnclvitllnp.  X,  Y,  and  Y1  value:;  ami 
their  associated  intervals  il.-tIotip,  with  interval  width  and  midpoint) 
lot  X values  across  each  'at  a ranp.c  can  he  t mind  in  the  attached 


<1 


listing.  The*  numerical  bohaviot  of  these  quantities  is  discussed 
below. 


Observa t ions : 

The  distribution  of  interval  widths  for  Y and  Y'  are  given  below 
for  BPA  and  BPF  versions  of  the  4-polnt  and  10-point  cases. 

1.  Distribution  of  interval  widths  for  4-polnt  case 
A)  Y intervals  for  BPA  version 


Interval  Width 

,„-n 


XXXXX  XXXXX  XX 

XXXXX  XXXXX  XXXXX  XXXXX  XX 


B)  Y intervals  for  BPF  version 


Interval  Width 
-11 


XXXXX  XXXXX  XXXX 

XXXXX  XXXXX  XXXXX  XXXXX  X 


C)  Y'  intervals  for  BPA  veision 


Interval  Width 
, -11 


XXXXX  XXXXX  XXXX 

XXXXX  XXXXX  XXXXX  XXXXX  XXX 


D)  Y'  Intervals  for  BPF  version 


Interval  Width 

,„-ll 


XXXXX  XXXXX  XXXXX  XXX 
XXXXX  XXXXX  XXXXX  XXXX 


II.  Distribution  of  interval  widths  for  10-point  case 

A)  Y intervals  for  BPA  version 

Interval  Width 
.„-7 


XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXX 
XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  X 
XXXXX  XXXXX  XX 


XXXXX  XXXX 


Observations:  (cent.) 


R)  Y intervals  for  BPF  version 


10-7 

XXXXX 

XXXXX 

XXXXX 

XXXXX 

10-8 

XXXXX 

XXXXX 

XXXXX 

XXXXX 

10~9 

XXXXX 

XXXXX 

XXXXX 

XXXXX 

io~10 

XXXXX 

X 

10'  11 

It)-12 

X 

0 

XXXXX 

XXXX 

XX  XX 

XXXXX  XXXXX  XXXXX  XXXXX  XXX 

I 


1. 


C)  Y'  intervals  for  RPA  vet  ; i.  n 

Q 

10  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXX 

10  "l>  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

in'0  XXXXX  XXXXX  XXXXX  XX 


01  Y'  intervals  for  RPF  v< rsion 
Interval  Width 

10  8 XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX  XXXXX 
10~9  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

io‘ 10  XXXXX  xxxx 


From  tho  tabulated  results  on  interval  width  distributions  Riven  in 
I and  II  above,  we  can  make  several  observations: 

a)  Comparing  4-point  > anil  Y'  widths  Y (for  both  RTA  tnd  RPF 
versions)  intervals  vat y insignificantly  in  v\llh  from,  co  re- 
ponding  Y'  Intervals  for  same  X values. 

b)  Comparing  10-point  Y and  v'  widths  same  as  for  a)  tbove. 

c)  Comparing  4 -point  Y ind  10-point  Y widths  widths  are  smaller 
for  4-point  case  than  for  10  point  case  (see  1A,  lit.  11  A,  lib 
distributions) . 

d)  Comparing  4-point  Y'  and  10-point  Y*  widths  - sane  as  for  c) 
above  (see  1C,  IP,  IlC,  TIP  distributions). 

e)  Comparing  4- point  Y for  RPA  and  RFF  versions  - widths  vary 
insignificantly  (see  1A,  IB) 

f)  Comparing  10  point  Y for  RPA  and  RVF  versions  - RPF  widths 
more  noticeably  wide:  than  in  4 point  case  (see  llA,  11R) 

g)  Comparing  4-pi  int  Y*  for  RPA  and  RTF  versions  - very  slight  in- 
crease in  width  of  RPF  intervals  over  RPA  intervals  (see  1C,  111) 


t 


Observations:  (cont.) 

h)  Comparing  10  point  Y*  for  BPA  and  BPK  vets  tons  more  noticeable 
increase  In  widths  for  HIF  over  UFA  than  In  4-point  case  (see 
IIC,  IID) 

2.  The  observed  Increases  In  Interval  widths  mentioned  In  comparisons 
c.d.f,  and  h of  observation  1 above  arc  reasonable  considering  that 
more  computations  are  performed  fi'r  the  10-point  case  than  for  the 
4-polnt  case,  thus  widening  some  of  the  computed  intervals. 

3.  The  widths  of  the  intervals  associated  with  Y and  Y’  (in  the  isngr  of 
10  ^ to  10  make  plotting  somewhat  imptactical.  Plots  (based  on 
3 significant  figures)  ot  the  uoulnterval  4-polnt  and  10-potnt  cases 
for  Y and  Y'  are  attached. 


Our  next  few 

observations  deal 

with  execution  time  of 

t he  Sl'l  INF  rout  lne 

MPA  time 

HPF  time 

HP  F t tme 

Noninterval  time 

Nonintorval  time 

HPA  t Into 

4-polnt  case  without 
overhead 

23.  15 

18.45 

.8  11 

4-polnt  case  with 
overhead 

h7 .2« 

5 b . .nl 

.817 

10-point  case  without 
overhead 

181  .71 

154.0/ 

. 804 

10-point  case  with 
overhead 

1 34 . 50 

108.05 

.803 

•< . From  the  tabulated  results  above,  we  see  that  the  HPF  interval  versions 
for  both  the  4-point  and  10-point  data  consistently  required  801  el 
the  execution  time  of  their  UFA  interval  counterparts  and  still  produced 
competitive  interval  widths. 

.*>.  The  noticeable  differences  between  the  ratios:  computed  with  function 

call  overhead  and  without  function  call  overhead  can  orobablv  be  attri 
buted  to  the  inability  of  the  SV.CON','  function  to  accurately  measu:  e such 
small  execution  times  (all  our  values  were  under  ' seconds!.  We  note 
that  the  l.itio  of  interval  to  noninterval  versions  Is  significantly  less 
for  the  4-polnt  case  than  for  the  10-point  case  and  there  are  fluctuations 
in  these  ratios  depending  upon  whether  o’  not  function  call  overhead  is 
Inc ludrd . 

b.  A more  consistent  execution  time  pattern  emerges  if  we  compart'  rat ros 
of  10-point  1 1'  4 point  CPU  values  for  each  version  tested. 


10  point  111’A  time  10  point  BPK  time  10  point  Non  interval  time 

f with  fen  call  ‘>.7  5.5  2.85 

ij  overhead 

i without  fen  5.7  5.5  7.00 

1 call  overhead 


J 

1 


| 


« 

■ 


•» 


Observations:  (cunt.) 


We  therefore  sir  that  both  interval  versions  (with  or  without  function  call 
overhead)  for  the  10-  point  case  require  approx inalely  5.6  times  the  execu 
t Ion  time  ns  their  corresponding  4 point  interval  version,  whereas  the  non 
interval  10-point  version  can  take  eithet  .'.S5  or  7 times  as  r.uch  execution 
tine  as  the  corresponding  4 point  non-intei*  il  time.  There  is  probably 
more  fluctuation  In  these  nonlnteival  it;..-  due  to  the  small  CI’U  times  in- 
volved and  the  inabllltv  to  mea  a rt  t>  prec  !•  ely. 

7.  The  results  from  the  WES  4-oitit  and  10  point  test  cases  do  not  seem 
to  strongly  justify  the  use  of  interval  analyst  h rt  ; the  midpoints 
of  the  Y and  Y'  Interval?  at  e approx  iriately  equal  to  tht  computed  non- 
interval  values.  Vosslbl  v ’..r:  '-'.Cunt  intervals  1 .•  might  re 

assure  the  user  of  a bound  above  and  below  the  true  ■ olution  for  some 
tihve  car  problem  but  wo  • est  ion  how  -a.  t ;k  an  ’ is,  situ  the 

given  data  was  to  no  nr  re  t w-  1 *. . >nl  ! mt  i igurt  and  it  is  doubtful 
that  such  figures,  in  pract  wo,  would  ho  exact.  frtotval  versions  of 
SPUNK  might,  however,  he  applicable  to  s L t cat  Ions  in  wh  h there  is  a 
spike  and/or  discontinuity  ir.  l u r.o.i  oared  a. it.-  in  wiere  a telero...ce 
(interval)  Is  given  to  ex  ri  i th<  uncertainty  Ln  h . ir  more 

data  points,  in  which  case  w * : g* t see  the  effect  of  si  ■ tolerances 
on  the  Intervals  associated  v Iti  tbo  use  o'  the  generated,  spline  liven 
for  such  a ease  with  one  o-  more  tolerarct  a lea. it  squat es  ipproaca 
would  probably  be  more  appropriate  t ' ? b-  rrobleiu  than  the  spi  ine  tech- 
nique. Further  testing  with  a wide  variety  ol  pract leal  data  seta,  of 
\ trying  sires  is  advised  before  deter-  -u Ing  the  ef £ ieaex  o!  me  interval 
version  of  Si’Ll  NK.  1 wou Id  also  sugcost  that  WES  co  le-  using  Alan 
Cline's  splines  under  tension  .uvpro  u (.see  "Scalar  - and  I'lanar  - 
Valued  Curve  Fitting  Using  Splines  Undet  elision",  Co-.  -,un  . at  Ion;  el  the 
Association  for  Commutin'’,  vac’- • net  v , Yo’ur.e  17,  Nw  ct  -,  Apt  . I lg  ■•) 
for  comparison  purposes  with  the  rent  ire  used  in  this  experiment. 


1 J . Problem:  To  sum  th«‘  first  l.''*  t« 
foiwurd  wltti  Interval, 
matte,  a tut  using  WPS  t« 

rms  ot  1 1 \) * * ( l l)  backward  and 
single  and  double  precision  arlth- 
st  data  tot  x values. 

Code  vised  and  references: 

A FORTRAN  routine  provided  bv  \,  t S 

Conditions  under  which  test  was 

omUn-  XeA: 

The  KTN  4 . h Kl^R  TRAN  compiler  was  used  with  no  opt  Imitation 
(l.e.  opt  “ O'.  Nonlnteival  (single  and  double  precision) 
and  Interval  (both  I'd  ' \ and  RPK)  eases  wore  run.  The  WRITF 
routine  was  used  with  vat  tables  declared  Interval. 

Ccmj>utat  tonal  Kesnl  t s 

Case  tor  x » 2. 

p 

Vo  war d S VU'.Y  at  ton 

Noninterval  (single  precision' 

1 . i'! 

Noninterval  (double  pi. vision) 

. OOOOOOOOOOOOOOOOOOOOOOOl'OOO  | 

Interval  (BPA) 

V l . oss  uioosusous.)  > .■ . 00 ''00000000 1 1 •< ' 

Interval  (RPV' 

i, . ‘N  M'-'i'i.w o-isss,  .' .0000000000000  i) 

Kackwat  d Summ.it  .mi 

Nonlnteival  (single  provision' 

1 anus )>)<) 

Nonlnterval  (double  precision) 

.’.OOOOOOOOOOOOOOOOOd  '00000000 

Interval  (BPA) 

( 1 . 0')0usS')0-o)0 o snu  t . O.'OOOOdiV'OOOOO' 

Interval  (BPKl 

^ l M i nio.MO'W'i'Mo , . Oi'OOOOOOOOOiV  " 

Case  lot  2 m -2. 

Vovwavvl  Sunm.it  .on 

Nonlnterval  (single  precision' 

. hboshbt'tiht'ht' s S 

Noni  ttiiv.il  tdouble  precision' 

.(■  '!•(>(>('  no  :.(>••(>  (>!'('('  oooooc-  t'pt't'  o’ 

Interval  (BPA' 

Interval  (BPV) 

in. 

’ 

BiiokwiU  v1  •vu.unat  Ion 

Nonlnterval  (single  preeiston' 

. C'C't't’t't' t't't' t'C'Ot'?'  » 

N'ninterval  (vloub'e  p>ee  Is  Ion) 

. (' t' t' f'f'f'f' f '(*(>('('(' t*  ’ 

Interval  (Id’A' 

Interval  (HIT' 

( . ‘w  , . * v> ) 

J 

L 

♦ t' 

-J-  i 

I'ompiit a t t on. 1 1 Results:  (.  c » * ; \ t . ^ 
C«80  fot  X .’.l 

N''ih  it  • i v a I v i n>*  I o pi  i t s i ..) 
Sonlnt  erval  ulowblv  \ 1 e<  1 . • <<n ' 
Interval  (Rl’A) 

1 nt  et  va  l (Pt'i  ' 

N v.il  v pr . 

Non  int. «■• v.\ l v,'.ou'-  !«■  | : eo  i ■ 
tntervil  (v'A' 

Int viva’  (RrP' 


Forward  Summation 
1 . 904090904090  r 7 
1.4  ' 04,l40‘»l  91  \ 1 

i -V  90"  •>1"ipiUV/,  l .409  "i.'-.,  ( i<) , t<S) 

' • 

Ha  ■ ' w.ml  Sia'iinat  1 on 
1 , 914  V4(.l')()4P909l 
1 . 4p4p-.«;0  >•  ' loo  1 17. 4 S'". U'.'SAt 
4040404(410909  1,1  .40  41'  >,<  i ' 40404.- 1 
1 .40-.  '4  <.1.1404040,  1.9090  U'4  ' • <4041  1 


to  e fot  x “ 0.1 

Non  ! nt  i"  \ a 1 (■  -ine  '.<-  precision'' 
Non  i i prei 

Intel,  at  (BPA) 

1 n t erva  l ( RPK ' 

Noni  , • \ \ ' (si ng , pi  ec  on) 
Non  • t i*i  v.i  1 (double  pie.  is  m\' 
Int  ere.  1 (SPA) 

Interval  (It  IT) 


veiw.ud  SoPBn.it  on 

o is-. 

. s ”41 " ' ssS  >S7p9p  l <■>'  t.S  < ' 1 1 S8 
( , 6 1 - • . - ■ < - 1 , , i ■ • ■ ' 5 ) 

. : * . . , ••  . , • .1 

Backward  Stirni.it  i i 
.t<’  7..19JS-.S  IS  7 0 *> 

. f<  ’ ■ . . <-.S  IS  ‘ 0 4 S N ■ 1 ’ ;.S  ■ " ’ 4». 

< . f ’ ’ • . 1 ■ ■ .S  >S  7 • , '•'■<  : " ' s' 

(.  . t-  ' 7.. ! 91  IS  ,l-  1 "I  • > ' " •<> 


Case  for  x *-  S.O 

Xonlnterv.il  (sin.  ’.<•  pree-sior.' 
Non  i nt  el  va '.  t iloub '.  < provision' 
l uteri  il  (SPA) 

Interval  (BIT) 

Non  interval  (single  pre.  is  ion) 
Noninteiv.il  (double  precision' 
Interval  (Rl’A) 

Interval  (BIT) 


Koi  ward  Su  - i.it  io-i 

1 . 4 999444000994 t 

1 .49944440. <<>9-144  >49  • '104  1444 

v 1 . 49O404<>9.)44 441  . . NOiV  ' !,  . 0 ' IS.'' 

i 1 . ■,444044449999  . . 1 . SOOO  ’Ot'i'O  \V  1 1' 

R.i.-kwa t Sun i.i  t Ion 

l . ,<)  }<'U<  090004  14 

1 NOi'OPOOOi  I’-'iV.'i'i'i'  ' '0  ' ' '0 

( l .49  94499409  44  >4.1.  xpoOdp  ' MOO.  V l ' 
l . •» 9 '9449449999,  i . M'OO.V  ' MOO.'l'l) 


Comput.it  tonal  Results:  (cent.) 


Case  tor  x “ -3.0 

Nonlnterv.il  (single  precision) 
Noninterval  (double  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BVF) 


Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 


Forward  Summation 
. 749999999999901 
. 749999999999901 

. 7A99 99999 9 99 99 999 99 9 9999999949 

( . 749999999999780, . 750000000000000) 
(.  749099999999901, . 7 50000000000 U 7 ) 

Backward  Summation 
.750000000000000 
. 7 500000000000000000000000000 
( . 749099999999996, . ' 50000000000004) 
(. 74  99  0999999  999 1 , . 7 50000000000004) 


Case  for  x * 5.0 

Noninterv.il  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interva’  (BPF) 


Forward  Summation 
1 . 24999999999990 
\ . 7 4 9 99  9 9 9 9 9 9 9 999  9 99 9 99 0 9 09  99 

( 1 . 999909900 00/ , 1 . 2 50000000 POOS 31 

(1 . 2 4 9 9 9 9 o o 9 o 9 9 9 7 , I . 7 50 00000000021) 

Backward  Summation 
l . 74999999999999 
1 . 25000000 0 0 0 000000000 0 0 0 0 0 0 0 

(1 . 27.999999999999  _ { . 7 '■OOOCOOOOO'OO  1 > 
(1 . 24990999999999 , i . 7 500000000000 1 ) 


Case  for  x » -5.0 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 


Forward  Summation 
.8333333333 3 3758 
.833  M33  33  '>33;n  i33  3 15  333  < 33  29 
( .8333  ) 1333  133105,  . 8 3 3 3 3 3 3 s 37. 3 i860) 
(.  833 3 VJ3  1 333  32 58  ,. 83  13  1 3 13 33 3 3.,p 3) 


Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 


Backward  Summat ion 
.8  3 1 3 3 3 33  33  3 3 379 


I 

i 


P 


. 8 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

( .8  < 3 3 3 1 3 3 3 33  3 3 3 7 , . 8 ' 3 3 3 3 3 3 3 3 3 3 3 3 b) 

( . 8 1 3 1 3 3 3 1 3 3 3 3 3 3 7 , . 8 3 3 1 3 3 S 3 3 3 3 3 3 3 b) 


Computational  Results:  (cont.) 
Case  for  x «*  7.0 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Case  for  x = -7.0 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Case  for  x =*  10.0 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 


Forward  Summation 
1.16666666666660 
1.166666666666666666666666666 
(1.16666666666660,1.66666666666751) 
(1.16666666666660,1.16666666666683) 

Backward  Summation 
1.66666666666666 
1 . 666666666666666666666666667 
(1.16666666666666,1.16666666666667) 
(1.16666666666666,1.16666666666667) 

Forward  Summation 
.874999999999936 
.8749999999999999999999999997 
(.874999999999769, .875000000000224) 
(.874999999999936, .875000000000057) 

Backward  Summation 
.875000000000000 
. 8750000000000000000000000000 
(.874999999999996, .875000000000004) 
(.874999999999996, .875000000000004) 

Forward  Summation 
1.11111111111106 
1.111111111111111111111111111 
(1.11111111111106,1.11111111111197) 
(1.11111111111106,1.11111111111126) 

Backward  Summation 
1.11111111111111 
1.111111111111111111111111111 
(1.11111111111111,1.11111111111111) 
(1.11111111111111,1.11111111111111) 
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II 


Computational  Results:  (cont.) 


Case  for  x « -10.0 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Non  interval  (double 
lnteiv.il  (BPA) 
Interval  (BFK) 

Case  for  x = 11.0 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Case  for  -11.0 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 


Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 


precision) 
prec i sion) 


prec ision) 
prec  is  ion'1 


precision 
precis  ion 


prec is  ion) 
prec is  ion! 


precision) 

precision! 


precis  ion) 
precis  ion) 


Forward  Summation 
.909090909090860 
.9040909090969090909090909089 
( . 90909090904068 3m . 40909090904 1137) 
( .904090409090860, .909090909090960) 

Backward  Summation 
.909090909090907 
.909.'  109090909090909090909091 
( . 909090909090907, .909090904090910) 
.909090909040907, .909090909090910) 


Forward  Summation 
1 . 09999999999996 
1 . 100000000000000000000000000 
1.09999999999996, 1 .10000000000087) 

( 1 . 09999999999996 , 1 . 1 000000000001 5) 

Backward  Summation 
1 . 00000000000000 
1 . 100000000000000000000000000 
(1,09999999999494, I . 10000000000000) 
(I.O9999999999999,  ’ . ’.00000  '.'000000) 


Forward  Summation 
. 9'.  66666666 .'661 c 
. 916666666666b  i't’66666'6666664 
( . 9 16666666 6 664 40 , . 0 166666 ''66667 1 4 1 
(.9166666b 66666 !8,. 916666 666666714) 

Backward  Summation 
. 916666666666664 
. 91 66666666066666666 66666666 7 
( . 9\  6b666666666t>4  ,. 916666666666668) 
( . 91666666666666* , . 916666666666668) 


r'0 


(com  . ) 


Ccmiput.lt  t Ollfl  1 Hc.-.u  Its: 

Case  for  x " 13.0 

Nonlntorv.il  (single  precision) 
Nonlutorv.il  (double  piool  si  on) 
Interval  (Bl’A) 

Inteival  (BPF) 

Nonlntorv.il  (single  precision) 
Nonlntervnl  (Joul >’ e precision) 
Interval  (lU'A) 

I nterva 1 (B" F) 

Cast  for  x “ -13.0 

Non  Interval  (single  precision) 
Nonlnterv.il  (dooMo  pr eo  Is  l on) 
Interval  (OPA) 

I nt  i rva  ' (BPF) 

Noulnteiv.il  (single  precision) 
Nonint«rval  (double  precision) 
Interval  (Bl’A) 
lntc*rvnl  (lU'D 


Forward  Summation 
1 . 0000000007692 > 

1 . 0000000007 69230769822485208 
(I  . 00000000076923 , 1 000001)00077014) 
(1  .000000000769/  3, 1 .00000000070925) 

Backward  Summation 

1 .0000000007092? 

1 . 000000000 769  .•  107698  .>  i/t  g ■,  >08 
( I .0000000007697  • ,1.0000000007692  1) 
( I .(IO(H)OUtKKl/b97  1,  1 .00000000076923) 

Forward  Summation 
. 99999992  10/687 
.99999999  7 30  7 (.9  7 16686  19048703 
(.99999999?  10746'),  .9999001)92  107920) 
( .990999992 107686, . 909900097  107696) 

Backward  Sumn1.1t  Ion 
.999009907  107686 
.999999997.107697  16686  19048  .’03 
( . 099999997  10  ,'689 . . 0000090. 7 10/69  1) 

( .9990 90  197  10/689 , .00000000  > (07 6°  i) 


O6serv.1t  ions  : 

1.  The  forward  mid  backward  double  precision  summations  usual.lv 
agree  or  differ  sllghtlv  (by  no  more  than  two  e.  clm.il  dlg'tsl 
w'th  the  backward  double  precision  rami  being  more  accurate. 

7.  Single  precision  summ.it lot"',  differ  bv  at  most  three  decimal 
digits  with  the  backward  sum.  being  more  accurate. 

1.  For  the  interval  cases,  Bl’A  bounds  tend  to  be  tig.htei  (or 
agree  with)  than  I1FF  bounds. 

4.  Bl’A  and  BPF  forward  sum  bounds  are  larger  than  UFA  and  BFF 
baekwaid  sum  bounds,  as  Is  to  be  expected  since  theie  is  more 
error  propagation  in  the  forward  summing  process. 

5.  With  respect  to  t ho  given  test  program  and  the  WES  data,  interval 
analysis  cannot  be  Justll  led  hete  on  our  CPC.  machine.  The 
double  precision  result  1-.  more  accurate  and  taster  than  interval 
analysis.  lhe  timing  i at  los  (til’ll  time)  of  double  precision  t o 
single  precision  summing  (both  backward  and  forward)  ranges 

from  ^ 1.25  to  ~ ?.'>,  wbore.i  the  timing,  ratios  of  interval  vet 
slims  (both  11FF  and  UFA  in  forward  and  backward  summing)  to 
single  precision  summing  ranged  f row  16  to  ^ 74. 


1 
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M.  Probl«:  To  subtract  I./2.*8  from  1.0.  repeatedly  square  the  dlf- 
In, mu-..,  and  »un  the  squate*.  Forward  and  backward  suras 
arc  computed  for  single  precision,  double  pi  op  is  ion  ,nJ 
in.nvs.  (Hr A ami  HPn  This  problem  was  j„  ined  by 

w s for  .hr  case  ol  1 V • . and  I have  also  rhos.n  to  , 
the  same  tost  with  1.  because  this  quant  i . v is  mo,  r 

appropriate  to  a CDC  machine  ami  roughly  analogous  to  \./2 
tor  the  I’nivae  1110.  the  mv  othn  write  up  fo,  details  of 
a I nt vac  111(1  case  run  on  a CDC  machine. 

Code  used  and  references; 

A FORTRAN  routine  hv  VFS  and  mo, lifted  to  include  double 
precision  case  and  use  tile  quantity  1.  7.sS  instead  ot  1 / ' ‘ " 

Conditions  undo,  which  test  was  conducted- 

Ihe  FTN  A.h  FORTRAN  compile,  was  used  with  no  optimisation 
O.e.  opt  - 0).  Noninte.val  (single  and  douMe  piecisfon) 
and  let  c.  val  (both  HV  A and  11  m case,:  we,  e The  WRITK 

toutine  was  used  with  variables  declared  lntcv.il.  The  st  in 
dard  (default)  option  ol  chopped  arithmetic  was  used.  The 
program  was  then  run  again,  using  the  civ  ROW,  (see 

observation  it  lor  cowpa, (son  of  rounded  and  truncated  results). 

1 at  i on.i  1 t osu  1 ( m : 

Non inlerva 1 (single  precision' 

Noninterval  (double  precision) 

Interval  (RVA) 

Interval  (BFi  ' 

Noni nt e, val  (single  precision' 

Noninterval  (double  precision' 

Interval  (HP A' 


lot  wa . d S um 

10. | o (,(,i  IS,;-  (.-poo  • ;oo  ts; 
t °,J  ■'  I ” f't'0  'M\  to  . oo.’  i ,j  • | j • i i 
( to.  oo.’  i o (,|,or  los  ...o. (’0011000000000) 

Ha,  It  ward  Sum 
Id  . It’-iOoo  i o -g 
<°.  d°.'|  o .'(,f,l  1 ft  1 1 - l.’poo  • -oo  p,  t\ 
(10  'JO.'  1 o '(.(,p  ' t O', > 1C  oo;  1 o •(,,,) , ,,  "I 


Nontnferval  (single  proofs  ton' 
Non i nt oi va 1 (double  precision' 
Interval  (HI'A) 

Interval  (Ht'V) 


lerms  oi  Sequent  o 
. OOOOOOOOOOOOOftO 
. 000, 10000, iqoooo.'gos.i  ■ 100(1 

( . 0000000,1000, too  1 > 00000000000000  It 

( . ooooo ooooooo oo t . . oooooooooooooo  ; t 


( 


Computational  Results:  (cont.) 

Nonintorval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPP) 

Noninterva1  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 

Noninterval  (single  precision) 
Nonintervnl  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 


Terms  of  Sequence 
.999999999999979 
.9999999999999857891 452847980 
( . 999999999999986 , . 999999999999986) 
( .999999999999986, . 999999999999993) 

.999999999999957 
.9999999999999715782905695963) 
(.999999999999972, .999999999999972) 
( . 999999999999972 , . 999999999999993) 
.999999999999915 
. 999999999999 94 3 L 3658 11391933 
(.999999999999943, .999999999999943) 
(.999999999999943, .999999999999993) 

.999999999999829 
.9999999999998863131622783899 
(.999999999999886,-999999999999886) 
(.999999999999886, .999990999999993) 

.999999999999659 
.9999999999997726263245567926 
(.999999999999773, .999999999999931) 
(.999 99 9999999773,. 99999999 99 9999 3) 

.999999999999931 
.9999999999995452526491136369 
(.999999999999545, .999999999999545) 
(.999999999999545, .99999999999Q993) 

.999999999998636 

. 9 999Q90 Q9 9 9 9090 505 29S 2 2 7480 5 

( . 999994999999091 , . 999999999999091) 
(.999999999999091, .999999999999993) 

.999999999997272 
. 9999999999981810105964557882 
( . 999999999998181 , . 999999999998181) 
( . 999999999998181 , . 999999999999993) 
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Computational  Results:  (cent.) 


10 


Noninterval  (single  precision) 
Noninterval  (double  precision) 
*)  Interval  (BPA) 
c Interval  (BPK) 


Terms  of  Sequence 
. 99*1999999994543 
.9990909909963620211929148851 
( . 999999999996162 , .999999999996362) 
(.999999999996362,. 999999999999993) 


Noninterval  (single  precision) 
j Noninterval  (double  precision) 
^ j Interval  (BPA) 

Interval  (BPF) 


.999999999989086 
. 09999999 09927 2404 2 38 584300 51 
(.999999999992229 , .999999909992724) 
(.999999999992724 .. 909999999999993) 


Nonlnterva'  (single  precision) 
Noninterval  (double  precision^ 
Interval  (BPA) 

Interval  (BPF) 


.999999999978172 
.99999909998544;  0847717389498 
( . 999999999985448 . . 999999999985448) 
( . 999999999985448  , . 909009999909993) 


Noninterval  (single  precision) 
Noninterval  (double  precision) 
\ Interval  (BPA) 
l Interval  (BPF) 


.999999999956344 
.9999999999708961695436896 
( . 9999999999 7 0896 , . 999909999970896) 
(.999999999970896, .999999999999993) 


j.  Noninterval  (single  precision) 

Noninterval  (double  precision) 

14  / 

Interval  (BPA) 

) Interval  (BPK) 

Noninterval  (single  precision) 

^ 1 Noninterval  (double  precision) 

] Interval  (BPA) 

's  Interval  (BPF) 

^ Nonintorval  (single  precision) 

, , Noninterval  (double  precision) 

16  J 

\ Interval  (BPA) 
i.  Interval  (BPF) 

Noninterval  (single  precision) 

Noninterval  (double  precision) 

\ Interval  (BPA) 
l Interval  (BPF) 


.999999999912689 
.999999999941/923390882263485 
( .999999909941792, .999999999941792) 
( , 999999999941792 , . 999999999999993) 
.909999999935377 
. 9999999999835846781798408288 

( . 9uo 99 999983 3 5 j . 999999999883585) 

( ,99999999988358 7 , .999999999990993) 

.999999999650754 
.9999999997671693563732341847 
( .99999999“767169, . 99999Q909767169) 
(.999999999767169, .999999999999993) 

.999999999301508 
.9999999995343387128006784779 
( .999999999534339, .999999999534339) 
( .909999999534339, .999999999909993) 
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Computational  Results:  (cont.) 


Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BI’A) 
Interval  (BPK) 

Noninterval  (single 
Noninterval  (double 
Interval  (BrA) 
Interval  (BPK) 

Noninterval  (single 
Noninterval  (double 
Interval  (lU’A) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPK) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPK) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPK) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 


prec  is  ion) 
prec i s ion) 


prec ision) 
prec is  ion) 

precision) 
prec is  ion) 

prec  i s ion'' 
precis i on) 


precision) 
prec is  ion) 

precision) 
precis! on) 

prec i s i on) 
prec ision) 


prec ision) 
precision) 


Terms  of  Sequence 
.999999998601016 
.999 9999 9 906867/425818 1973 902 
( . 999999999068677, . 99999999906867  7) 

( . 99999999906867  7 , .999999999999991) 

.999999997206032 
.9999999981373548525037565177 
( .999999998137355, .999999998137355) 
(.999999998137355, .999999Q99999993) 

.999999994412065 
.9999999962747097084769599808 
(.999999996274710, .999994996274710) 
(.999999996274170, . 99999099999999 J) 

.9999999^3824129 
.9999999925494194308317077177 
( .'199999992549419,  .999999992549419) 
(.999999992549419, .949999999999993) 

.999999977648258 
59999998  50988389171 74  5662  531 
( . 99999998 50988 19 , . 999999985098839) 

( . 99999968 50988 39, . 660999999999993) 

. 99qqq99552965l 6 
.0999999701976780563937341224 
( . 999999970197678, .990969070197678) 

( .9999999970197678, .999909099999993) 

.999999910593033 
.999999940 395357000965861 47 57 

( . 499999940395335, . 999999940305359) 

( . 099999940395355,  .99969999999999  '•) 
.999699821186073 

.0666668807907 l 7 55464 5 18966. ’8 

( .999999880790714, . 699969880760721 ) 

( .690966880760714, .999999990999993) 


Computational  Results:  (cant.) 


26 


28 


29 


10 


31 


32 


33  J 


Non  In rerval  (single 
Mon interval  (double 
Interval  (Bl’A) 
Interval  (.BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (RVA) 
Interval  (HP!1) 

Noninterval  (single 
Non  interval  (double 
Interval  (BFA) 
Interval  (BPV) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (RFC) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (HPK) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (RFF) 

Noninterval  (single 
Nonlnterv.il  (double 
Interval  (BPA) 
Interval  ( BP V ) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  IRI'K) 


Terms  of  Sequence 

precision)  .999999642372174 

precision)  .999999761  .4814493201434011219 

(.999999761481439, .999199761481446) 
(.999999761481439, .999999999999993) 

precision'  .999999284744472 

precision'  .9999995231629554836921105270 

( .999999523162931 , .999199523162970) 
( . 99999952  11  (<2931  , . 99999999  9999993) 

precision)  .999918469481451 

precision)  .9999990463261383409512441013 

( .999999046 326085, . 999999046326167) 

1.999999046326084. . 999999999999993) 

precision)  .9999971 30980948 

precision)  .9999980926531861757368998850 

( . 19999809264 3076. . 99999S0n2653243) 

( . 9 9 9 1 1 80 1 • 7. 6 5 30 7 6 , .999999999999903) 

prec is  ion)  . 9999942 77970078 

precision'  . 9 9 9 o 9e 1 8 5 3 1 00 1 0 3 2 33  4 2 00  5 3 3 8 1 

( .999096184 1 09787 ,. 999 9 96 IS. S3 101 24) 

1 . 9999961 85 30978 7 . .  909099990999093) 

precisian)  .99999855.4972894 

precision)  .94Q912  3706  345  7250640134997  73 

( .999992 3706391 21 , . 999992370634800) 
1 . 999992 3 7 Op  34 121 . .911199909999993) 

precision)  . 9999771 12076741 

precision)  .9999847413273522296289345357 

(.999984741326447, .999984741327808) 
( . 99998474 1 376447 , .999999999999903) 

precision)  .999944724677  360 

prec  is  ion)  .9999694828874  3 1 5.402296847379 

( .999969482884721 , . 999969482888442) 
( .999969487884721 , . 9Q999919991999 3) 
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Computational  Results:  (cont . ) 


Non  Interval  (single 
Noninterval  (double 
Interval  (RPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Non Interval  (double 
Interval  (RPA) 
Interval  (BPF) 

Noninterval  islngle 
Non  Interval  (double 
Interval  (RPA) 
Interval  (RPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (RPA) 
Interval  (RPF) 

Noninterval  (single 
Noninterval  (double 
lnter\ al  (RPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 


precision) 
prec Is  ion) 


precis  Ion) 
precision) 


prec is  ton) 
pr  eels  Ion) 


precision'' 

precision) 


precis  1 on) 
prec  islon'1 


precision) 
precis  ion'' 


precision) 

precision) 


Terms  of  Sequence 
.999908451450100 
.9999 189b670b357253871 381 9093 
(,99993896670273b, .999918966708179) 
(.9999 'S9b6702736, .999999999999993) 

.999816911281316 
.O9Q877O371377774406244388169 
(.999877937130535,. 999877937141402) 
( .990g/793 713053 5, .9999909990949^3) 

.990633856084149 
.0997558891748972152123903751 
( .O99'5588oi60411  , . 90o ; SS889182183) 
(.999755889160411, . 999000009990993) 

.999267846229664 
.9995118379398893627871643624 
( .999511837910923,  .990511  837950458) 
( .099511837910923, .999999900999993) 

.098536228508470 
.9990239141819756570357081981 
(.909023914124070, .999023914211101) 
( . 0090239141 24070, . 999900009999093) 

.90707459964392 0 
. 908048781 1074754623221848 234 
( . 9°8 048 780991 777 , .998048 78116Rb~2) 
( . 9O80487809O1777 , . qooqo 900900 090 3) 

.994157757255081 
.09610136047701174b'  1276869595 
(.906101360239170, .996101360086285) 
( .oog 101 360239170, .999000009000903) 


Observations: 


1.  Backward  and  forward  double  precision  sums  (noninterval) 
agree  to  1<*  decimal  digits,  as  do  backward  and  forward 
single  precision  sums  (noninterval). 

2.  Both  single  precision  sums  (backward  and  forward)  agree 
to  3 decimal  digits  with  double  precision  result.  (This 
compares  with  9 decimal  digit  agreement  for  WES  case  using 
\.H.  instead  of  1./2 .‘■•S.) 

3.  The  lower  bound  for  both  interval  versions  of  the  sum 
agree  to  3 decimal  digits  with  the  single  precision  ver- 
sions for  backward  sums.  (This  compares  with  15  decimal 
digit  agreement  for  WES  case  using,  1./2.27  instead  of 
1./2.48.) 


4.  The  lower  bound  for  both  interval  versions  of  each  term 
in  the  sequence  agrees  to  a steadily  decreasing  number  of 
decimal  digits  (13  digits  to  2 digits)  with  the  single 
precision  results  for  ail  40  terms.  (This  compares  with 
15  decimal  digit  agreement  for  WES  case  using  l./2.^  in- 
stead of  l./2.^*8.) 

5 The  BPF  intervals  for  each  term  are  equal  to  or  larger  than 
the  corresponding  BPA  intervals. 


6.  The  BPA  intervals  for  the  terms  of  the  sequence  are  de- 
generated through  term  23.  (This  compares  with  only  two 
degenerated  intervals  for  WES  case  using  1./2.27.) 

7.  The  number  of  decimal  digits  in  agreement  for  upper  and 
lower  interval  bounds  on  the  terms  of  the  sequence  drifts 
downward  as  follows:  15  digits  (up  to  term  23),  14  digits 
(term  24),  13  digits  (term  25),  12  digits  (term  28),  10  di- 
bits (term  30),  9 igits  (term  38),  8 digits  (term.  39). 

(This  compares  with  a decrease  from  15  digits  to  0 digits 
for  WES  c;>5?  with  l./2.-?i,  but  for  that  case  decrease  started 
by  term  4,  whereas  here  decrease  doesn't  start  until  term  24 
and  never  goes  below  8 digits.) 


8.  As  each  term  in  the  sequence  is  computed,  the  upper  and 

lower  interval  bounds  are  menotonicly  decreasing,  i.e.,  the 
sequence  of  lower  interval  bounds  is  steadily  decreasing 
from  terms  1 to  40  and  the  sequence  of  upper  interval  bounds 
is  also  steadily  decreasing  from  terms  1 to  40.  The  resulting 
intervals  for  the  terns  of  the  sequence  are  drifting  to  the 
left,  however,  unlike  the  WES  case  with  1./2.27,  the  inter- 
vals for  l./2.^®  terms  do  not  collapse  to  near  zero;  this 
is  because  the  difference  between  1 and  l/2^8  is  much  closer 
to  1 than  the  difference  between  1 and  1/2^';  hence  powers  in 
the  former  case  are  larger  than  in  the  latter  case. 
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Observat ions : 


(cone . ) 


9.  The  single  precision  resultant  of  the  sura  Is  not  contained 
In  any  of  the  computed  interval;  tot  the  sum.  (Recall  for 
WES  case  with  1./2.*-'  it  was  the  double  precision  resultant 
of  the  sum  rounded  or  chopped  to  single  precision  that  was 
not  contained  in  any  of  the  computed  intervals  tor  the  sum. 

10.  The  number  of  decimal  digits  in  agreement  in  the  corres- 
ponding upper  and  lower  interval  bounds  for  the  sum  is 

13  and  15,  respectively.  (This  compares  with  13  [for  both 
bounds)  for  the  WES  case  for  1./.’.-'.)  For  each  interval 
for  the  backward  and  forward  suns  there  is  n 10  decimal  digit 
agreement  between  upper  and  lower  hounds.  (This  compares 
with  8 decimal  digit  agreement  In  WES  case  for  1.  1 

11.  When  the  program  was  run  with  the  CDC  rounded  option  we  get 
the  results  given  in  the  following  table  (truncated  results 
and  interval  results  are  listed  for  comparison). 


Noninterval  (single  precision  truncated' 
Noninterval  (single  precision  rounded) 
Interval  (UFA) 

Interval  (RFF) 


Forward  Sum 
39 . 98830s 0991 07b 
40 . 0000000000000 

(39.9‘>."lu  t.t>0;i9.S,  }0.0QJin?pt, I.U71) 
( 39 . 99  2 1 9 t.607 19 .5 ,40 . 0000000000000) 


Noninterval  (single  precision  truncated) 
Noninterval  (single  precision  rounded) 
Interval  (PEA) 

Interval  (RFF) 


Backward  Sum 
39. OSS  1040991078 
40.0000000000000 

(39.9021 97P607  JOS  ^ . qq;  i q 70*6 14 1 62' 
( 39 . 99 ; i o ' (-S071O *, , 40 . 0000000000000' 


We  note  that  the  noninterval  single  precision  rounded  and 
truncated  lesults  (for  both  backward  and  forward  sums)  are 
net  In  the  BPA  intervals  for  their  respective  sums  tin  cen- 
ttast  to  WES  case  for  1./2.  > although  the  rounded  single 

precis’ on  results  are  contained  in  the  RFF  intervals  for  sums 
The  single  pr  vision  rounded  terms,  unlike  their  single  pr e 
c is  ion  trtineated  counterparts,  are  all  exact  In  equal,  namely 
1 . 0000000''0000000;  thus  for  the  rounded  single  precision  ease 
all  the  differences  are  erroneous  (rounded  to  1)  and  thus 
the  rounded  sum  of  40  is  also  inaccurate. 

12.  I strongly  suspect  for  this  test  with  l./2.^®  on  out  CDC 

machine,  that  the  best  result  is  the  UFA  interval  associated 
with  the  backward  sum.  The  overall  results  for  this  tost 
are  probably  slgnl  • ie.rntlv  worse  than  the  results  for  the  WF.S 
case  with  1./2.*',  because  in  out  case  with  1./2.48  we  are 
pressing  the  roundoff  limits  of  out  machine  for  the  powers  of 
the  difference  between  l 1/2^*^  ti.e.  1 l./2^)  and  Its  powers 
ate  much  more  inaccurate  that  1 1 l./2*-0  and  its  powers,  rspc 
ciallv  since  our  machine  has  a 48-hit  binary  mantissa). 


?7 

Problem  To  subtract  1./2.  from  1.0,  repeatedly  square  the  difference, 
and  sum  the  squares.  Forward  and  backward  sums  are  computed 
tor  single  precision,  double  precision,  and  interval  (BPA 
and  RIF)  cases.  This  problem  was  defined  bv  WF.S. 

Code  used  and  references: 

A FORTRAN  routine  provided  bv  WES  and  modified  to  Include 
double  precision  case. 

Conditions  under  which  test  was  conducted: 

The  FTN  4.6  FORTRAN  compiler  was  used  with  no  optii  izatlon 
(i.e.  opr  « 0).  Non  ntevva'  (single  and  double  precision) 
and  Interval  (both  BPA  and  BPV)  cases  were  run.  The  WRITE 
routine  w.is  used  wit!:  variables  declared  interval.  The  stan- 
dard (default'  option  o*  chopped  arithmet*  • was  used.  The 
program  was  then  on  avail  using  the  • oun;  option.  See 

Observation  II  to  comparison  of  rounded  and  truncated  re- 
sults. 


Computational  result . 

Noninterval  (single  precision' 
Noninterval  (double  precision' 
Interval  (RPA) 

Interval  (RPV) 

Nor  interval  (single  precision' 
Non Interval  (double  precis  ion) 
Interval  (8PA) 

Interval  (RTF) 

Non interval  (single  pro. is • >n) 
Noninterval  (double  prec  sion' 
Interval  tBPA) 

Interval  iBPF' 

Noninterval  fr ingle  precision) 
Non  interval  (double  prod*  ton) 
Interval  (BPA) 

Interval  (PPF) 


Forward  Rum 
2S.bt>7  ’s;’N  .PS  ?6  ’ 
r''.ps?-'s:p.’l  ",  R.'’Rp0p"  .07184 
( . ' ' *»3  ■>..  pS  167,.' 1 s . bb  ~ ' r<2t  • ' N -»31 ' 

t,2f>.  -b72:>254bS.tb  ' , 7' . 6b  7 2 5 f.°l  55  S' ) 

Backward  Sum 
2 5 . pb ^ 7 s 7 S4 1>S3  1 
2s.  t-b  S7b21 ' lc." 800b 7 7407 184 
i, 7 s . bb  ' 7 S .'  ‘-..nS  ' T 1 , . bf'  ’ 3 S 7b  S41S) 

< 2S  . b - 152  whSlM.  .'s.bb'7 54191.5555) 

Terms  of  Sequence 
. Qoocoougspogg  to 
. ooo  ' ooog  .('os.R  lS8Cb‘  **73.« ' ' s 
(.ooqqqqogsoosB.Io,  .0.'  00th  ' S7RR83  •' 
^ _ 9000000s  6048839 , .oqooqco  • ,p  )■ 

, 909004970197078 
. 9099909  7019  b '834  '•.03074350 
( . 99994  99  1)1 9 7b  7S  , . 90999997019  7 b’S' 

^ _ 9990999 7019767 S, . 9990000 ; 014  7 b86) 


Computational  results:  (cont . ) 


NonlntervaL  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  IBPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPP) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 


precision^ 

precision) 


prec Is  ion) 
precision) 

precision) 

precision) 


precision) 

precision) 


precision) 

precision) 


precision) 

precision) 


precision) 

precision) 


precision) 
precis  ion) 


Terms  of  Sequence 
.999999940395355 
.9999999403953565568769913153 
( . 999999940395355, . 999999940395359) 
(.999999940395355, .999999940395377) 

.999999880790714 
.9999998807907166664675026124 
(.999999880790714, .999999880790721) 
(.999999880790714, .999999880790760) 

.999999761581439 
.9999997615814475437882381192 
(.999999761581439, .999999761581456) 
(.9999997 61 581439,. 999999761581538) 

.994999523162931 
.9999995231629519309826315538 
(.999999523162931, .999999523162970) 
(.999999523162931,. 999999523 163137) 

.999999046326085 
.9999990463261312355356742820 
( .999999046326085, .999999046326167) 
( . 999999046326085, . 999999046326504) 

.999998092653076 
.9999980926531719649193127447 
(.999998092653076, .999998092653243) 
(.999998092653076,. 999998092653243) 

.999996185309767 
.9999961853099819017610409730 
(.999996185309787, .999996185310124) 
( .999996185309787, .999996185310124) 

.999992370634121 
.9999923706345156634562602886 
(.9999923706341 21,. 999992370634800) 
(.999992370634 121,. 999992370637 525) 


Computational  Results:  (cone.) 


Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (P"A) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BFF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Nonlnterva’  (single 
Noninterv.il  (double 
Interval  (UFA) 
Interval  (BPF) 

Nonlnterva I (single 
Non  interval  (double 
Interval  (BPA) 
Interval  (B°F) 


Terms  of  Sequence 

precision)  .9999847* I32b447 

precision)  .Q999S4  7*13  272  38b4460ol063M7 

( .949984741  <23>**/. . 999984 74 1 1 2 7808) 
(.999984741326447, .99998*74133 1261) 

precision)  .99996948288572k 

precision)  .9999694828873041836933934933 

(. 90946048288 V21 , . 999460482888  . .2) 
(.9999e°*S2f'S5??l , .ooqo.v 1482899353) 

precision)  . 999938°66702736 

precision)  .0999389667059025345959422444 

( . 0494 389667027  36,  99°o  3.c'466708174' 
<,  .9494  <896670.  736,  .99uoagot,(,~ 30003) 

precision)  .990877937130535 

prec  t sionl  . 4998; 79  37136869057  5 795898 196 

(.0998  7 79 3 ;i 30535, .9498  7 ;o3?14i  4 20) 
(.94987'c  17130535, . 009877937 1 S5072) 

precision)  . 909' 85889160411 

precision?  .0907358891730786711264978204 

(.  9997559  591604 1 1 ,.  999  7 5588'*  182 183) 
( . coo  :3599c : 60411 , . 9Q0  7SS8  4 162 IS 3) 

precision)  .999511837910923 

precision)  . 9Q95U8379362531624670871601 

( .99951183  74;  03  3 . .999511937954.'.  -S' 
( .994  <118179  :0923, .99951  ;r :8129013) 


Noninterval  (single 
Noninterval  (double 
Interval  (BrA'' 
Interval  (BFF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 


precision' 
pr<v  i s ion) 


precis  ion) 
precision) 


.909023914124070 
. 090(37 )0 1 4 ' 747068065073846316 
( .99907 3914124070, .9990239'  .211101) 
(. 999023914124070 , .999023914560045) 

. 998048  ■'80Qol  777 
.998049781092951951  7004195°  •'  3 
(.998048780991777, . 998048781862881) 
(. 998048780991777,. 9980487P18b288l) 


Computational  results:  (eont.) 


Noninterval  (single  precision) 
/ Noninterval  (double  precision) 

^ Interval  (BPA) 

i 

^ Interval  (BPF) 

Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 


Terms  of  Sequence 
.996101369239170 
.9961013694411271244005132349 
(.996101369239170,. 996101369586285) 
(.996101369239170, .996101370977986) 

.992217937800149 
.9922179382024888262313622151 
(.992217937800149, .992217938491674) 
(.992217937800149, .992217938491674) 


21 


Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 


.984496436092378 
.9844964368907979353036467808 
(.984496436092378,. 984496437464667) 
(.984496436092378, .984496442966631) 


Noninterval  (single  precision) 

J Noninterval  (double  precision) 

/ 

) Interval  (BPA) 

Interval  (BPF) 


.969233232678594 
.9692332342506768817987186284 
(.969233232678594, .969233235380624) 
(.969233232678594, .969233246213953) 


r 

23  / 

] 


Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 


.939414059328594 
.9394130623760274857324544856 
(.939413059328594, .939413064566395) 
(.939413059328594,-939413085566443) 


Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (EPA) 

Interval  (BPF) 


.882496896037107 
.8824969017627061076775694322 
(.882496896037107, .882496905878028) 
(.382496896037107, .882496945333472) 


f Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

^ Interval  (BPF) 


.778800771515126 
.778800/186207753543801752689 
(.778800771515216,-778800788884293) 
(.778800771515126, .778800858522914) 


26 


Noninterval  (single  precision) 
Noninterval  (double  precision) 
Interval  (BPA) 

Interval  (BPF) 


.606530641712553 
.6065306574531306230190265813 
(.606530641712553,-606530668766798) 
(.606530641712553,. 606530777236031) 


63 


Computational  results:  (cont.) 


C Noninterval  (single  precision) 

J Noninterval  (double  precision) 

| Interval  (BPA) 

^ Interval  (BPF) 

/■  Noninterval  (single  precision) 

Noninterv.il  (double  precision) 

' Interval  (BPA) 
v.  Interval  (BPF) 

f Noninterval  (single  precision) 

Noninterval  (double  precision) 

^ Interval  (BPA) 

V.  Interval  (BPF) 

Noninterval  (single  precision) 

Noninterval  (double  precision, 

' Interval  (BPA) 
g Interval  (BPF) 

Noninterval  (single  precision) 

Noninterval  (double  precision) 

Interval  (BPA) 
interval  (BPF) 

Noninterval  (single  precision) 

Noninterval  (double  precision' 

Interval  (BPA) 

> Interval  (BPF) 


Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 

Noninterval  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPF) 


precision) 

precision) 


precision) 

precision) 


Terms  of  Sequence 
. 367879419336241 
.3678794 3843052687 18 784263932 
( . 367879419336241 , . 367879452154700) 

(.  367879419 J 36241,. 367879583734547) 

.135335276171169 
.1353352812199598161617013234 
(.135335267 1711 69, .135335291317642) 
(.135335267171169, .135335388128705) 

.0183156345 

.01831  ii<3S3i  288560810628670080 
(.0183.15654  5402917,  .0183156410760311) 
t.‘  183156345402917,. 0183156672799473) 

.000335462668613526 
.000335467607907381 4645397657185 
( .00033546246861)526, .000335462703025999) 

( . 000335462468613526 , . 000335463667909734) 

.112535067848281  X 10*6 
.1125351613040215530993168484  X 10*6 
(.112535067848281  X 10-6, . 11253522847613  X 10  6) 
(.117  )506  8281  X 1 0“6, . 1 1 25358724 37453  X 1 

.126641414956172  X 10-13 
.12664167.52972214994343244580  X 10~13 
(.126641414956172  X 10*13, .126641776481763  x 10* 3 ) 
(.126641414956172  X IP"3"3 , .126643225965123  X 10~13) 

.160380479821014  X 10~27 
.1603810125792185243493947999  X 10*27 
(.160380479821014  X 10*27, .160381395504569  x 10*?  ) 
(.160380479821014  X 10*2 7 ,. 160385066328534  X 10*27) 

.257218983076186  X 10~55 
.2572206919595345054490025814  X 10~55 
(1.257218983076186  X lO*35 ,. 257221920239931  X 10*  5) 
(.257218983076186  X lO*55, .257233696615933  X 10*55) 


64 


Comput.it  tonal  results:  (emit.) 


)S 


36 


17 


3 H 


l'l 


Nonintvrval  (single 
NonlnLerval  (double 
Interval  (BPA) 
Interval  (BPK) 

Noninterval  (single 
Non  Interval  (double 
Interval  (!U’A) 
Interval  (BPK) 

Noninterval  (single 
Nenintetv.il  (double 
Interva  (BPA' 
Interval  (BPK) 

Nonirtevv.il  (single 
Noninterval  (double 
Interval  (BPA) 
Interval  (BPK) 

Honin’ erval  (single 
Nonintcrv.il  (double 
Interval  (BPA) 
Ir’erv.ll  (BPK) 

Non  In’,  ervn ' (s  ■ ng  ’ e 
Non  1 u* erva t (double 
Interval  (BPA) 
Interval  (BPK) 


prcc Is f on) 
prec is  ion ) 


Terms  of  Sequence 
. 66 1 6 1 6057  54  74  72  X ]0~',S 
. 66 162484  i / 2049 1 3952 73  5181 01 2 X 10 


-55 


prec  is  tun) 
prec i s Ion ) 


prei  i s i on) 
pree l s f on) 


pi  ec i s l on) 
| * ’’  < ■ . ' : s i -Ml  1 


pi  ec  1 ::  ion) 

pi  e<  >n ) 


over  i s ie •:  1 
pr ec i s f on) 


(.  1 s;.":89830761S6  X 10’  .257221920219931  X 10  S5 
(.257218983076186  X 10”  25723369661 6933  X 10~V’ 

.437/  35800988497  X lO-^ 
.43774743382816466208049256??  X 10-*** 

■ .43  77  l '800°88497  X 10  **‘,.4  1775579  7164S0  X 1 ' 
;.<*!/  *35800988497  X 1.0  * * *' , .437835967  ’ ! »ost.*  X 

0. 

0. 

_ ?oa 

(0. , . 1 >65/5653123/01  X 10  ) 

-i’ll 

(0.,  .156575653.175701  X 10  * 1 


0. 
0. 
(0. 
(0.  . 

0. 

0. 

(o.. 

(0., 

0. 

0. 

(0.. 


-■>()< 


..156575653125701  X 10 
.’.5657  56 53 175 7 11 ' X 1 0~  *’ ° 1 


. to  a 

. 1 5657565317 5701  X !0  * ) 


. I 56' 


5s  531 


- 20  l 
U X 10  ) 


.1565756^31  7 •.':»!  X !•) 


-?o  t 


(0.  , . 1 56 5 '56 5- 1 75  -01  x 1 a '3 


Ba t'kw.i re!  and  f oi  ar.!  ■'  1'  ! pi  ec ts  fo*t  'uuns  (non ' nterva ! ^ agree 
exact  I y whi’nM  Vick  we  • 1 and  ■ oivar-.’.  s(  eg *.  e prec  i - on  ‘>uws  y non 
interval)  agr  •«.  to  l.l  dccfrvil  dip.  t$  ot* t of  ' S. 

2.  Both  n ingle  precision  sum.  (S;i  -kward  and  forward)  agree  to 
9 decimal  dig fts  with  double  precision  result. 

1.  The  lower  hound  for  both  Interval  version.-;  of  the  sum  agree 

exactly  with  the  single  precision  version*,  for  backward  and  for 
ward  sums. 

4.  The  lower  hound  for  both  Interval  versions  of  each  term  lit  the 

sequence  agrees  with  th«  single  precision  result  for  all  40  terms. 


t»f3 


Observations  (coni.) 


Non  ♦ tit  «• 
Non  Into 
l nr  * rva 
Into?  v*i 


S.  The  RPF  intervals  for  both  sums  uni  tot  each  term  are 
equal  or  latter  than  the  or  spending  UFA  intervals. 


The  number  ol  decimal  digits  (not  including  leading  zeros) 
in  agreement  for  uppei  uni  lower  interval  bounds  on  the  terms 
ol  the  sequence  dr  lit-  downward  a-  follows:  IN  dibits 
fter  ms  l and  7 ) , l l digits  item  4 t l dibits  (term  7), 

1 2 dibits  (term  ’),  11  1 1 g i t s (term  ll),  1 0 digits  (trim  14), 

n digits  (term  17),  8 digits  (term  70),  7 digit*  l term  2D, 
b digits  (term  10),  S digits  (term  H),  4 digits  (term  14)  , 

0 digit:;  (term  17)  . 


7.  As  each  term  in  f'\e  *«c<iu«mk  e i>  computed,  the  upper  and  lower 
interval  hosad  a »a*  mono i on  i.e  decree  mu-,  l.e..  the  sequence 
of  lower  ■ i hounds  Is  stead  Iv  decreasing  1 1 om  ferns 

* to  40  arid  : h«  sequence  of  vrpper  interval  bounds  r also 
ste.-idilv  d»s  • ag  1 i in;  ?»  rs  1 to  »0.  The  t*-\  If  ng  inter- 
va  ! r;  t or  • « t • • . . of  the  • ‘ quer.  • i * e d ’ ft  g t O the  1 e<  t 

appr  oaeh  o-g  no.  v i na  l l v , • he  ntetv;i’  e • «*  r t i . ' ’ l v collapse 
from  term  14  onward. 


g The  double  pr *i  is  ion  rrult  .mt  o*  the  sum  it  r Minded  or 

chop  ped  ' u;l  e pi  •*.  on  s eon?  a nod  n im  ol  the  com- 
pute 1 interval:;  ! or  the  sum. 

0.  The  number  oi  decimal  digi'  i n • .«  s*u(  i r the  cor  responding 
upper  and  lover  nterv.  ’ bound*  to*  * he  sum  ;s  13.  For  each 
interval  ! r he  backward  rnd  or  ward  sums,  there  Is  an  8 
dec  In  «•  dlg.it  a* reenenf  et  v<  en  up"«u  an.  lower  hounds. 


10.  Accord  ng  to  the  comments  n this  Wl  S listing,  the  results 
lot  the  suht  r i lion  on  he  thrive  l MO  we*  •'  correct  using 
interv.D  but  net  exarf  without  interv.nl  date  t o truncation 
or  rot  ' r M < * ’oaf  i eg  -'oic  wonder  1 'u'eival  I't’sul  t was 

;>,'t  t ■"  , not  entirely  because  of  true,  at  ion  problem.  but  be— 
causi  • nf  rv.'!  imp  l orient  a * po:.s  v w *1  double  precision 

on  the  1 0 wh  ch  l be l e\e  has  a wider  exponent  range  than 

s f ng  l e prec  (si.'  1 and  r hu  » ov  hi  • < g t • a.  k o!  sma  1 l er  qu  rut  i 
tier  . rat  ha  than  set  them  equal  to  ; -mo  as  hr  Single  preci- 
sion f ! on  ; ’ 'V  po'of  a t t f turret  r • . On  *N-  1 ’ 1 P version  of  this 
test  were  double  •'  • .»n  ”»n  interval  arithmetic  results 

f or  t h ’ • 'ir  ii‘,  . t ei  « . H t ’ ■ rrt  e*  v r ’ re*  u 1 1 s 


U. 


When  f he  program  w» 
tire  result*  given  in 
and  in!  <m  v.j!  ii-.uh  s 


run  with  the  CPf  rounded  option  we  get 
ie  o*  owing  f.dh  (trunafed  re  nit 
are  M * ed  for  compnr  Ison): 


Forward  Sum 

t'ol  ( ••  ingle  pre  is  ion  truncated)  ' > . bb  / 7 V' N4b8  lb7 

i\  a l (single  precision  rounded)  '>.H-  7N7N89SO/8 

l (ft  PA)  ( 7 ** . (>b  7 '2  S 2 >4  bS  167,  7 % ,6b  7 ? N 7 b 7 N S4  1 1 ) 

\ (\\yy\  ( . bb  7 S4b8n  M , S . bb  7 7 S 11  ^ l‘>  S s 7 ) 


Observations:  (cont.) 


Backward  Sura 
25.6672525468373 
25.667252J895078 


Noninterval  (single  precision  truncated) 

Noninterval  (single  precision  rounded) 

Interval  (BPA)  (25.6672525468373,25.6672526755418) 

Interval  (BPF)  (25. 66725254b8373 , 25 . 6672531915555) 

We  note  that  the  noninterval  single  precision  rounded  and 
truncated  results  (for  both  backward  and  forward  suras) 
are  in  botli  BPA  and  BPF  intervals  for  their  respective  sums. 
The  single  precision  rounded  terras,  although  differing 
slightlv  from  the  corresponding  truncated  versions,  exhibit 
the  same  monotonic  decrease  to  zero  that  the  unrounded  terms 
ind icated . 

12.  Although  not  provided  wit*'  the  exact  results  from  the  llnivac 
1110,  I would  strongly  suspect  that  on  our  CDC  machine  the 
best  result  is  the  RPA  interval  associated  with  the  backward 


15.  Problem:  To  use  Kalian's  problem  Involving  roundoff  error  propa- 
gation to  check  out  consistency  of  compiler  generated 
arithmetic  rules  versus  execution  time  arithmetic  rules 
and  to  determine  the  effects  ot  the  application  of  interval 
analysis  to  this,  problem. 

Code  used  and  reference;; : 

We  used  as  a basis  for  our  test  the  FORTRAN  code  given  by 
F.  Oort  and  C.  Mole.,  S I ON CM  Newsletter,  Volume  8,  Number  2, 
April  1973,  pp.  2 4-26.  This  code  Is  as  follows: 


Program  1 

Program  2 

i^Nr  - 1.0 

TW0  - 2.0 

THREE  - 3.0 

FIVE  ■»  5.0 

H - 1. 0/2.6 

11  - iANE/TVii? 

X *>  2.0/3.0-H 

\ - TW0/THREE-H 

V - 3.0/5.0-11 

Y - THREE/' FIVE  -H 

K = t X* \ 1 X) -H 


y r: 

(Y+Y+Y+Y+Y) -H 

Q - 

?.0*F/F. 

PRINT, Q 

Note 

• 1: 

The  code  In  the  middle 

belongs  to  both  Program 

is  1 and  2. 

The  correct  value*;  (»> r 

E and  F are  both  Eero, 

which  means 

that  Q should  be  indeterminate:  however,  on  most  (if  not 

all)  machines,  the  v.ilu 

o ‘ ot  Q is  t in  I to  ! For 

example . 

results  for  0 on  the  CPC  6600  and  CDC  .600  ns 

reported  by 

Dorr  and  Molet  anil  on  t 

he  CPC  Cyber  70.  model 

72  as  report 

hv  Cinsberg  are  us  toll 

ows : 

Program  1 

Program 

CDC 

6600 

(truncated  arithmetic) 

-6.0 

+3.0 

cnc 

6600 

(rounded  arithmetic) 

-6.0 

-6.0 

CPC 

7600 

(truncated  arithmetic' 

6.0 

+ .3.0 

CPC 

7600 

(rounded  arithmetic) 

-6.0 

+4.0 

CDC 

Cyber 

70.  model  72 

-6.0 

+ .3.0 

(trunc.it  ed  ar  1 1 hmet  i c ' 

CDC 

Cyber 

70,  model  72 

+4.0 

-6.0 

(rounded  arithmetic) 

Not. 

' 2: 

The  difficulty  with  thi 

; problem  is  that  the  e 

valuation  of 

arithmetic  subexpression  cm  he  in  a different  order  and/or 
use  dll (erect  rounding  procedures  if  done  at  translation 
time  than  the  hardware  would  perform  at  execution. 


t>8 





No to  2: 


(cont . ) 


Furthermore,  this  inconsistency  is  compounded  by  the  sub- 
tract (ve  cancellation  and  roundoff  error  propagation  intro- 
duced throughout  the  problem  (especially  in  expressions  for 
E and  F) . 

Note  3:  Ideally,  we  would  like  the  application  of  interval  analysis 

to  ttiis  and  similar  problems  to  automatically  warn  the  naive 
user  that  reported  finite  values  of  computed  results  (c.g. 

Q in  tliis  test)  can  he  in  serious  error.  We  would  like 
interval  analysis  to  warn  us  in  botli  Programs  1 and  2. 

Conditions  under  which  test  was  conducted: 

The  FiN  4.6  FORTRAN  compiler  was  used  with  no  optimization 
(i.e.  opt  “ 01.  Noninterval  and  interval  (both  BPA  and  BPF) 
cases  were  run  with  RClUNl)  arithmetic  and  truncated  arith- 
metic. IN'IVR  was  used  to  output  interval  variables.  Also 
for  nonlntervnl  tests,  both  rounded  and  truncated  arithmetic 
options  were  used  under  each  compiler  option  (i.e.  opt  " 

0,1, 2. 2, VO) . 

Observat ions : 

1.  Tlic  results  are  as  follows  on  our  COC  1'vhor  72  witli  opt  - 0. 


Values  of  Q 


Program  1 

Program 

n 

a 

Von in ter va 1 ( t run)  -6.0 

+3.0 

Noninterval (rdd . ) +4.0 

-6.0 

I nterva  1. (BP A)  (-6.0, -6.0) 

(-.1265014083  10‘3?3,. 

.1265014083 

103?3 

Interval (BPF)  (-6.0, -6.0) 

(-.126^014083  1 O'333.. 

. 1 265014083 

103?  ' 

Program  1 is  a test  of  compile  time  arithmetic  behavior 
Program  2 is  a test  er  execution  time  arithmetic  behavior 


2.  From  the  results  given  in  Observation  1 above,  we  see  that 

BP A and  BP F interval  results  agree  for  each  program;  however, 
both  interval  versions  were  deceived  in  Program  1 because  the 
interval  routines  use  ?./’!.  and  3./5.  as  arguments  which  per- 
mits tlic  compiler  to  do  the  arithmetic  here,  and  unfortunately, 
the  compiler  does  not  evaluate  these  express  Ions  and/or  use  the 
same  rounding  rules  as  are  used  at  execution  time.  Thus, 
compile-time  arithmetic  deceived  the  interval  analysis  into 
indicating  the  true  result  for  Q is  -6,  when  in  reality  Q 
should  be  0./0. 

3.  The  intervals  for  ()  in  Program  2 were  not  deceived  because  no 
compile  time  arithmetic  was  involved  in  computing  tin-  intervals; 
such  is  the  case  when  interval  analysis  is  applied  to  arithmetic 


6‘  > 


f 

! 


(cent  . ' 


Obsei  vai  i ons  • 


expression  where  all  value  are  hidden  behind  variable  names 
and  in'  explicit  arlthinetl  wlrch  mav  be  done  in  several  dif- 
ferent wavs  Is  present.  thus,  for  Program  the  interval 
bounds  for  0 Implv  that  the  value  ot  q can  be  between  -00  and 
+ 00.  As  a mattei  el  tact,  l’ti'i  i an  ' Is  baited  bv  interval 
package  bv  warning  ot  pcs-  I bio  division  by  -ero.  Thus,  in 
Program  1',  tnvolv.ng  no  c.  -pile  time  aritbmet  . •,  the  interval 
bounds  proper Iv  warn  tin  u -or  that  the  result  lei  O is  pio 
bably  in  set  ions  error. 

Prom  (bis  test  am!  the  co-nont.-  in  Observation  we  see  that 
when  ! ,u  i'd  with  i - oinpti ; at  1 ana  l problem  in  which  we  intend  te 
use  Interval  analvs  s,  we  should  re  wr i t « (or  avoid'  all  ex 
pliclt  float  ing-polnt  •>'  • thmet  ic  express 'ons  (nuke  them  •’  ■ 
plleit,  i.e.  hi-'e  all  l.i.a!  tng  point  number ■ hei-  ml  variable 
names)  which  max  be  evaluated  bv  the  comp. lot  in  a d.lfietent 
manner  that  the  li.trdwaio  e e at  exeeutien  time.  This  pee 
cautien  will  decree-  e the  uiet's  vh.iaoc  of  h-'in-  deceived  -v 
the  resulting  interval  bound.-..  It  w ' 1 not,  howex-et  , c..u.tnfeo 
that  no  doce-t  • on  will  ut  . 'Yc nli.nl ti es  In  leu  iw.no 
hax.v'i  mil  attribntes  oi  interval  -ia v-.  . vlw.tvs  allow  room 
f or  the  introduction  i>!  tlecopl  ve  >ounds 

b . Tin'  nonlnt  i-'v.i  1 results  reported  ;n  0be-.vat;  nv  1 tor  Ptoc.rans 
I and  2 * tol  on  tent  ( ■ both  >und«s  and 
vorvion-O  be  -,-ise  of  in  otiststen  ies  hc'vi-i'it  c --p-  let  t ime 
and  execution  tee  artfhmo  ic  w'licb  .m-  fin  thet  magnified  ''V 
roundof  eti-':  pi  opa  -a  t leu  induced  l'\  subtractive  - aneel  1 at  1 on 
and  ton -existence  ot  exact  t ioat ing  point  representation  for 
some  of  the  quantities  involved  ye  g.  1.  S . ' . 

b.  The  noninterval  version  ot  Progi am a ’ and  ' were  run  en  out 
Cx’bei  72  va-.iier  opt  iove  l' , l nc  b.-' h -oimdod  and.  .‘hopped 

arithmetic)  a-  wo  l ’ a-  un.l  t pm  ave  • ••'  • lion:  with 

rounded  and  sopped  arithmot  »e ' The  re  n’  s tot  ’ . V,  and  q 
are  eixoti  below  ' ot  each  . as.  • 


Nonlnti’rva ' Results  on  Cvber 


OPT-O. 

Trim 


OPT-O. 

Rdd 


OPT-1 , 
Trim 


’’ru-hiim  1 


v . xss?;  M Tp’SiSi'd  s'  y ’ p 

F -.106581410 3b  .01 Xll'  1 
q-b . 00000000000000 

l . 3 5 52 ’1  3678300  'OX  10  ' 
k . : , dw ; 3S7KHOO X 10  ! 
Q 4.00000000000000 
l-  . 355271367880050  XlO  ’ 

f .lobbsi-'iioib  ioisy  io  ’ 
q t< . oooooooooooooo 


Pt  oc.v  am  ’ 


. 710'.-.:.  't-oioov  io 
. . On  '81  '*10  36-e'l  X It'* 
3.O0OOOOOOOOO0OI' 

. f-x;  i tt.  sst'ox  'X  io 

. 10658U10  KwOl  h x lo“ 

b. 00000000000000 


• 7l0r' 


i - bOiooX  io  ’ ’ 


. I 0b 481410  It'-'iOl  •'  )(  10 
1 . 00OO00000O0000 


Program  ' 

. ■’  Hlb  4 d • 3S760100  Xlu 
. 10b’>bl4103b401sK  li  " 1 ' 

3.00000000000000 

. ihs.’.’i  in  'dsoosoK  io 

. Itlt-hSls  '.  Olt'AOl  'X  '< 

6.00000000000000 

- .71054.'  • 1S/('0100  Vli 
. 10('SKlsl03bs01  'X’  io 
' . 00000000000000 


Noninterval  Results  on  Cyber  72  (cont.) 


0PT=1, 

Rdd 


OPT*=2 , 
Trun 


0PT=2 , 
Rdd 


0PT=2, 

UO.Trun 


0PT=2, 

UO.Trun 


OPT=2 , 
UO.Rdd 


Program  1 
. 355271367880050 
. 710542735760100 
4.00000000000000 
"E. 355271367880050  10 

F-  .106581410364015 
„Q-6 . 03000000000000 
.355271367880050 
.710542735760100 
4.00000000000000 
.355271367880050 
.710542735760100 
4.00000000000000 
.355271367880050 
• .1065814010364015 
•6 . 00000000000000 
.355271367880050 
.710542735760100 
4.00000000000000 


10 

10 


-14 

-14 


-14 

IQ'13 


10 

10 


•14 

-14 


10 

10 


r14 

-14 


10 


-14 

-13 


10 

10' 


-14 

-14 


Program  2 
.355271367880050 

- .106581410364015 
-6.00000000000000 

- .710542735760100 

- .106581410364015 
3.00000000000000 

.355271367880050 

- .106581410364015 
-6.00000000000000 

- .710542735760100 

- .106581410364015 
-6.00000000000000 

- .710542735760100 

- .106581410364015 
3.00000000000000 

.355271367880050 

- .106581410364015 
-6.00000000000000 


10 

10 


-14 

-13 


10 

10 


-14 

-13 


10 


-13 


10 

10 


-14 

-13 


10 

10 


-14 

-13 


Program  3 
.355271367880050 

- .106581410364015 
-6.00000000000000 

- .710542735760100 

- .106581410364015 
3 . 00000000000000 

.355271367880050 

- .106581410364015 
-6.00000000000000 

.355271267880050 

- .106581410364015 
-6.00000000000000 

- .710542735760100 

- .106531410364015 
3 . 00000000000000 

.355271367880050 

- .106581410364015 
-6.00000000000000 


10 

10 


-14 

-13 


10 

10' 


-14 

-13 


10 

10' 


•14 

-13 


10 

10' 


-14 

-13 


10 

10' 


-14 

-13 


10 

10 


-14 

-13 


A Program  3 was  added  which  is  identical  to  Program  2, 
except  the  constants  1,2, 3, 5 are  read  in.  The  interval  version*- 
of  this  test  were  not  run  under  round  and  truncated  options 
(and  under  other  opt  levels)  because  they  would  have  no  numeri 
cal  effect  on  the  interval  versions,  since  interval  arithmetic 
is  written  in  assembly  language  and  the  round  or  truncate  op- 
tion available  with  the  CDC  FIN  F0RTRAN  compiler  only  is  ap- 
plied to  FORTRAN'  explicit  arithmetic  expressions,  not  to  assembly 
language  code  or  library  routines. 


7.  Observing  the  results  given  in  Observation  6,  we  note  that 

a)  All  truncated  results  are  consistent  for  each  program 

under  all  options  and  all  rounded  results  are  consistent 
for  each  program  under  all  options,  but  all  Program  1 
and  corresponding  Program  2 results  differ  for  all  cases, 
and  are  not  consistent  for  rounding  and  truncating  under 
the  same  option. 
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7.  (cont.) 

b)  Corresponding  results  of  Programs  2 and  3 are  consistent 
under  all  options. 

c)  The  results  for  program  1 and  corresponding  Program  2 
under  truncated  and  rounded  options  is  identical  pattern 
for  all  opt  = 0,1,2,2UO 

d)  Although  Programs  1,2,3  look  equivalent,  we  see  that  we 
obtain  values  of  -6,4  and  3 for  Q when  the  exact  result 
should  be  0/0! ! ! 

8.  Although  the  CT)C  6600,7600,  and  Cvber  72  all  have  the  same 
f loating-pcir.t  arithmetic  format,  both  h’otc  1 results  above 
and  Observation  6 results  indicate  differences  in  values  of 
Q on  the  same  machine  as  well  as  on  different  machines. 
(Particularly  obv  ms  is  how  the  Cyber  72  result  on  Program  1 
differs  'ron  Progr.nr.  1 results  on  the  other  machines). 
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16.  Problem:  To  tost  the  performance  of  ■- 1 1: . cos,  and  tan  routines 
with  iutctv.il  analysis  as  specified  by  WFS. 

Code  used  and  references. 

A l 0 Rl KAN  program  provided  by  Wi S which  output  interval 
results  via  use  of  INTRAP  routine  In  interval  analvsis 
pac kage . 


Conditions  under  which  rest  was  conducted: 

The  FTN'  4.6  FORTRAN  compiler  was  used  with  no  optimization 
(i.o.  opt  0).  Both  SPA  .and  BPF  interval  versions  vtere 
run.  Interval  variab’es  wee  output  via  use  of  INTRAP 
routine  in  interval  package. 

Observat ions : 

1.  Both  BPA  and  BPF  versions  gave  identieal  results  (fot  both 
decimal  and  correspond  i -ag  octal  values)  f or  corre-  tid  ing 
interval  bounds. 

7 . All  iuterva’  endpoints  were  checked  foi  accuracv.  All 

cot  sect  1 1 i > v.t  \ :■■.  and  corresponding  arguments  are  indi- 
cated hv  double  arrows  on  the  output  provided  in  rite  at- 
tached listing. 

1.  As  suggested  hv  a comment  in  the  VI'S  to.  t deck,  we  checked 
endpoints  aganst  adjacent  computat  i ons . The  arrows  on 
our  listing  indicate  the  correspondences . There  was  no 
turthei  suec  t'ic  indication  ot  the  pattern  of  correspon- 
dences. Presttmab'. v t ie  pattern  observed  is  correct.  (We  . 
know  the  co- rospond.ences  indcatoJ  botv.-en  arguments  and 
endpoints  via  double  arrows  are  correct  because  these  were 
ver i f ied . ) 

4.  A few  modification:  wort  made  in  the  WVS  test  deck  before 
the  prog:  ;rt  was  run.  In  the  VS  vers. on  of  subroutine 
NTRAP , the  arrays  MON , NAMES.  TVPA.  TYB.  TYFk  in  COMMON 
block  IN;  it'.  vero  rack  specified  as  a single  array  of  di- 
ntens  i on  watch  disagrees  w.tb  flit*  dimension  ot  these  ar- 
rays a specified  Sv  V.  he  in  his  interval  package  . Thus. 

I changed  the  dimension  of  each  of  the  tforeraent  toned  erravs 
to  40  to  be  in  agreemen"  with  Yohe's  packag. . Also,  this 
test  deck  requires  the  sc  of  subvout  no  ' \ TVNS  v,1or  the 
INION  operator^  ;n  the  Yohe  package.  bt  : we  were  not  pro- 
video  with.  : tic.t  routine  in  the  version  given  to  us.  1 
have  written  a version  of  INl'VNN  (see  attached,  listing)  which 
satisfies  \ohe  s del  ini  t ion  o*  the  behav'oi  of  the  Oil  ON 
operator . 

h.  Since  BPA  and  BPF  interval  versions  produced  the  same  results, 
1 have  at t ached  only  the  BPA  listing.  It  shows  all  the  ar- 
gument . and  interval  endpoints  (in  decimal  and  octal)  pro- 
duced bv  INI  RAP  and  also  provides  a copy  ot  t ho  test  program 
ns  run  on  our  Cyber  72. 
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17.  Problem.  To  use  the  non-interest  and  interval  versions  of  the  WKS  FKT 
program  with  sin  input  data  provided  by  WKS. 

(.ode  used  and  references: 

KK1  pr°!ram  ts  wrltu'"  1,1  FORTRAN  and  based  on  the 

£ ‘FrS'M1:  jts  sr 


' 3*1 '■  1 *22®  onuo:  wh ich  t est  was  condt 


v.  V:'' v"  ; 1 '■■■- used  **«  » oPt u on  < . 

, ' : ' ‘ « va  anc 

. , . ' ' n*  ' oT'nev  were  or >ta  -n  bv  using  [he  S'  U‘V' 

function  to  “leisure  t ’io  C”,  * i*ne  > , 

i ' r . ‘ up  ludliK  t i mi*  { t’r 

• ■ ' ■ M • U rvai  vai  , . 

I ' ’ ■ ■ VLl  to  rep]  eth,  I 

GT  and  r,  respectiv.  . N 

_ **'  0 optput  valu<  nterval  variables.  I v,  \"  ■ . e. 

s ' 

' - w ; ■ i oVZl]  ; 2 ■ To? 

; ■ . . “ . -- 

l>  Use  the  right-  ..  

!<l>  • : ’ >««  «d  to  b< 

absolute  value  (see  Observation  ..i 


1!u‘  inTut  test  data  was  provided  br  a WKS  vvr  ar.vt.r  m 

2 ' ' ?6  ‘ ““  ts  of  an  array  A ..  •.  the  1 , el* 

1 ■ the  array  is  represented  bv 

•VI ( I ) ' S!\  (/ MEGA  * S(l" 


«TA  “ ATANdd.  .o.  ) * 800. 

A ' - n.OAT(l-)  ) S DT 

].  (F’,i\\T(\>  * PK' 

OK  = . S 
N - *00  fi 

F rout  -re  hould  -e  found  in 
'■  " ■ < ' ■ ores  accept  o,  , , 1 * ei 

■ » • to  X(  200)  where  in  the  FT  dt  v« eg, 

x(n  - float  u-n  j>  py 


Comp  tit  at  ion  a ] rose. I t s • 

(A)  Timings  for  the  KFT  Test 


Computation  results:  (cont.) 


Noninterval  version  ^2.8  seconds  of  CPU  time 
BPA  interval  version  v.'5l  seconds  of  CPU  t inn 
BPF  interval  version  '".'Is  seconds  of  CPU  time 


(B)  Computer  values  of  \(l)  corresponding  to  X(100) 
Noninterval 

. 49999  99999  9^988 


Note: 


BPA  interval 


(.-W999  99999  97368,  .50000 


00000  02125'' 


midpoint  of  interval  = .49991  99999  l9Q0(i 

-10 

length  ot  interval  = . 2685851541 1 7318  * 10 

All  the  values  el  Y(I)  for  1 - 1 to  2049  art1  given  for  non- 
interval  version,  all  values  of  Y(l)  other  than  r . one 
stated  I , - on  thi  >rdei  >1  ' . L0"  . 

For  the  BPA  version,  all  Y values  other  than  the  one  in  B are 
on  the  order  of  ' and  for  5PF  version  on  the  order  el 

lO-ld. 


Observat  ions : 

1.  Two  o!  the  statements  in  the  Wl'S  modification  for  line  1160  in 
the  ITT  cod*  are  : ".correct  as  given,  namely 


kT(L2>  = \ 

and 

WP(Ll)  = Y 


In  tlie  section  of  code  immediately  preceding  line  .1160  V is 
being  computed  to  the  necessary  powers  and  uler's  formula  is 
being  used.  c/e/.  e ■ co  ■ z -*■  1 s n . "lie  mistake  made 
in  the  UTS  modification  was  to  put  the  real  part  of  this  rela- 
tionship in  r'ae  imaginary  oca'  ion  and  the  imag  nary  part  in 
the  real  loeat  ion  in  the  array  WP.  Thus  the  correct  mod  if  .ca- 
tion should  be 


WP(Ll)  * X 
wp(i.:n  = y 

If  these  changes  are  lot  made,  t lie  code  produces  several  non- 
zero values  for  Y(I)  ind  not  at  the  expected  place  (i.e., 
corresponding  t-  X(200)) 

From  the  computational  results,  we  see  that  the  noninterval  and 
interval  versions  produced  values  close  to  the  exact  result  ot 
rhe  noninterval  version  is  slightly  less  than  the  exact 
value  o'  . and  both  interval  versions  included  the  solution 
tas  well  as  the  noninterval  resultant'.  As  in  our  other  tests 
with  the  interval  package,  t h<  BPF  version  produces  a slightly 
wider  nterval.  We  also  note  that  the  midpoint  of  the  BPA  re- 
sults is  les;  than  t it  exact  value  but  the  BPF  midpoint  is 
si  lghtl\  more  than  the  exact  value.  For  both  the  BPA  and  Bt'F 


Observations:  (com.) 


interval  versions,  the  intervals  for  the  other  values  of  1 
were  small  nonnegative  neighborhoods  of  zero;  all  ot  these  in- 
tervals for  BP  A and  BFF  versions  were  bounded  by  exact 1\  zero 
on  the  lett  and  the  lengths  i these  intervals  were  on  the 
order  of  10  H and  It  ' >r  BPA  and  BPF  versions,  respectively. 

3.  Examining  the  timing  results  lor  this  test  we  note  that 

BPF  interval  yer - ion  CPI1  t ime 
’A  interva  verslo  PU  t inn 

BPA  int  erva  was  on  CPU  time 
Noninterva.  v.rsion  CPU  time 

BPF  Interva vers  en  epu  t i»e 
Xi  ti Interval  version  CPU  t ime 

!hus  we  see  that  both  interval  versions  require  substantially 
more  CPU  time  than  ; e non interval  version  and  that  the  BPF 
interval  version  requires  only  S3  of  ’he  CPU  time  needed  by  the 
BCA  interval  vorsi  . and  provides  an  interval  result  or.lv 
slightly  larger  than  the  BPA  interval. 

Fart  of  the  .VS  mod  If  icat  ion  tor  1 ine  310  in  the  rrT  driver 
program  included 

Y ( I ) = SQRT(AFS  * ArR  + .VF!  * AFI) 

When  t u interva . vor-iicc.s  at  tht*  V FT  r,'s*  wer.  run  with  this 
oxoression  for  ' ’.'3,  he  urograms  halted  with  . . lag  'Stic  in- 

dicating that  the  program  wa.s  attempting  to  calculate  the  square 
tv'ot  of  a negative  u uitv.  c understand  why  this  '.a  V . 
recall  that  t ne  interval  .;c*  in  it  ton  for  the  mult  " c it  an  '• 
tve  iiit  erva  Is  ; 


■n in  (a*e. 

a*d. 

b*c , 

b*d'' . 

Mx  (a*c. 

i*d. 

b*c . 

b*d  ' 

when  both  intervals  are  the 

same , 

i.e.. 

a.b  = e.d  then 

a ,b  * a,b  - m.  i n (a 

, a*b , 

b“) , 

max 

. a*b.  b~) 

It  turns  out  that  for  most  of  the  YU)  cases,  the  le‘t  end- 
points  of  t ;;e  intervals  for  \FR  and  Ar  l are  negative  (relatively 
sr.ia  . . on  the  order  ot  : 0 - : ! v and  the  right  endpoints  are  pcs  i - 
tlve  and  sc  that  rult  Iptving  such  interval:-  by  ; ’-.er.se ives  does 
not  produ.-c  positive  Interval  . i.e. , Lf  a 0 and  b 0 
then 

i.b.  * . n,b  ~ ,a*b  i*b,  max  (a'  .b*) 


which  would  contain  negative  numbers  since  a*b  ■ 0 and  thus 
would  make  the  interval  square  root  routine  halt  because  it 

7b 


.83 


89.  t- 


'b . 4 
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recognizes  that  the  interval  is  not  entirely  nonnegative.  To 
correct  thit;  situation  we  rewrite  ill)  as 

Y (I ) - SQRT  AHS  (AFR  * APR  + AF  1 + AF  D ] 

which  results  in  tin  interval  version  of  the  square  root  routine 
working  with  an  entire.lv  nonnegative  interval.  This  situation 
is  a direct  result  of  the  rules  tor  interval  arithmetic,  i.e., 
an  interval  multiplied  b\  itself  does  not  necessarily  produce  a 
nonnegative  resultant  interval  even  though  for  noninterval 
arithmetic  on  rea  numbers,  the  oroduct  of  any  real  number  with 
itself  is  nonnog.”  . vc . 

■.  In  the  v.‘S  Inpu:  'or  the  FFT,  4096  values  of  t ie  SIN  routine 
were  used  with  arguments  between  0 and  ' 800  . S'nce  in- 

crea singly  ' argot  an  ; les  were  used  and  errors  in  the  SIN  rout  trie 
tend  to  got  bigger  as  the  angles  increase  beyond  2-  . further 
oxaniaat  n:  ' • SIN  results  was  conducted.  he  values  of 

•Yu  ,e.S  SIN  emits  wer>  cor  arec  with  . ie  double  precision  SIN 
of  the  double  • recistin  expansion  o'  the  original  single  preci- 
sion arguments  and  with  the  resultants  chopped  to  single  preci- 
with  ie  WES  IS  . >■  t ion,  the 

KVS  S'N  ar ie  it  expressi  ms  were  ca  ’ ulated  in.  doub  e precision 

am.  Y eti  used  with  doubl  proc  talon  81 N with  the  rest  ts  chopped 
to  single  precis' on  for  oompar : son  with  t'ie  «T8  SIN  values.  For 
hot  h sets  c*  co" oar isons  wt’o  the  WF.S  SI'.'  v.  : u s,  it  was  noted 
that  there  was  a si  ight  gradual  loss  of  accuracy  in  the  SIN 
values  .!  *,  he  angles  increased  in  sice;  t ere  was  never  more 

than  • oi  i it  Loss,  ...  erroi  novel  um 

or.  the  rder  o:  10“  l1’  for  the  abovenent  ioned  oompar  •'  sons . 0? 

with  >ns , we  are  tat  i it  the 

COt"  KTN  version  of  DSIN  is  more  accural  e tor  larger  angles 
(when  the  double  precision  resultant  is  chopped  ;o  single  pre- 
cision'1 than  he  sin  routine  is.  An  examination  o!  the  In- 
tervals i s l ’c  it  c ’ w i tit  tile  elements  in  the  urra\  whic  i 
contains  Ye  i . uitatit-  o*‘  t he  WKS  SIN  input  s also  reflects  this 
loss  ot  dig'ts,  ....  the  lengths  of  the  intervals  associated 
w th<  elemen  iri  y S tncreai  * slightly  as  the  sin 

angles  et  larger;  ar  the  RFA  version  the  interval  lengths  are 
initial  0 ar  on  the  order  of  0“  ' am  eecrease  to  the  order 
of  10-  -t  or  10“ 10  as  the  arguments  of  the  SIN  approach  800- 
and  for  the  liPK  version  • e ieterv.ils  are  initially  on  the  order 
of  10'1-5  and  decrease  to  the  order  of  10-l°  as  the  angle  ap- 
proaches 800i  . 

n.  More  am1,  varied  eases  would  >vo  to  be  run  boy  one  the  NTS  test 

>ef ore  wo  could  make  more  cor  his  ve  assessments  ot  the  behavior 
of  the  interval  versions  et  tb  WKS  FFT  pri'gram  with  respect  to 
exact  results  and/oi  nonlnt erval  runs.  \lso  it  would  be  in- 
formative tv’  perform  Interval  ami  non  interval  tests  on  ot  her 
tiigli  quality  FFT  routines  such  as  that  provided  by  the  IMS1 
I ibrarv  or  Subrv’utine  VYVKT  by  Nv'rman  Hi  oncer  f vlescr  ibevl  in 
"Three  Fortran  Programs  that  ’’ertorm  the  t'ooley-Tukov  Fourier 
Transform,"  Techn  . al  \ctv  l^c MIT  Lincoln  laboratory, 
l.exingtou,  Massachusetts,  .tu.lv  78,  1967  talso  available  from 
NtlS  under  AD  657  0191'. 
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APPENDIX  IS:  INTRlJDUCTION  TO  ERROR  HOUNIH NO 
TECHNIQUES  AND  INTERVAL  ANALYSIS 


In  this  appendix  the  reader  is  Intri  need  to  interval  analysis  as  an 
error-bounding  technique.  We  discus  the  motivation  for  error-bounding 
methods  and  give  some  examples.  We  t’-en  move  on  to  a discussion  of 
engineering  applications  and  the  advantage:;  of  using  such  techniques. 
Finally,  we  present  some  interval  mad  vs  Is  definitions,  theoros,  and 
elementary  properties.  Limitations  theoretical  and  practical)  of  in- 
terval analysis  are  described  ii  Section  of  this  report.  7>  > reader 
desiring  addil  tal  nformacion  • nterval  analysis  is  urged  t<  onsult 
Moore  (1966)  and/or  any  of  the  references  given  in  the  bibliography  of 
this  appendix. 

Most  scientific  and  engineering  problems  require  solutions  of  an 
extensive  range  of  mathematical  equations  - algebraic,  transcendental,  in- 
tegral, differential,  etc.  Quite  often,  an  approximate  or  numerical  solu- 
tion is  needed  because  an  analytic  or  exact  solution  is  not  known  or  Is 
computationally  awkward  to  work  with.  Such  approximate  methods  introduce 
their  own  inherent  roundoff  or  truncation  errors.  Thus,  the  scientist  or 
engineer  who  utilizes  these  techniques  frequently  would  like  to  know  the 
relationship  between  the  approximate  and  exact  solutions,  especially  when 
dealing,  with  math  models  used  in  designing  various  physical  objects  such  a- 
airplanes,  electronic  equipment,  or  nuclear  power  plants. 

Unfortunately,  many  numerical  solutions  give  little  or  no  indication 
as  to  whether  they  are  above,  below,  or  oscillate  about  the  exact  solution. 
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For  example.  Table  1 and  Table  2 from  Ralston  (19t>5)  show  some  output 
from  approximation  techniques  such  as  Hamming's  Method,  Milne's  Method, 
and  the  Hermlte  Method  when  applied  to  the  following  two  initial  value 
problems : 

(a)  y'  - -y  y(0}  - 1 

(b)  y’“  --y-  y(0)  - 0. 

1+tan’ v 

Scattered  throughout  these  tables  (underlined  entries'  are  changes  in  sign 
between  two  adjacent  entries  in  the  error  columns;  th  -so  places  represent 
positions  at  which  the  approximate  solution  esc  llates  about  the  exact  so- 
lution. Figure  1 from  Laptdus  and  Seinfeld  (1°'1>  reveals  that  the  fifth 
order  Nvstrom  approximat ion  f o'  problem  (a)  above  is  Initially  equal  to  the 
exact  solution  at  x«T',  goes  below  the  exact  solution  for  0 * x < -o  u.d 
then  rises  above  the  exact  solution  for  x > *5.  This  type  of  alterr.  i t ing 
behavior  is  quite  commonly  exhibited  by  numerical  solutions. 

A more  pronounced  example  of  approximation  fluctuation,  with  respect 
to  an  exact  solution,  is  exhibited  by  Greenspan  (1''67).  Three  possible  ap- 
proximate solutions  to  an  initial  value  problem  Involving  the  l;ne.i;  ,-ed 

(i)  ' , vl' 

Kmden  differential  equation  are  v and  v'  In  Figure  2 and  v in  Figure 
The  exact  solution  is  shown  by  the  sol  d line  curve  in  Figures  and  ?.  it 
is  apparent  from  these  graphs  that  ail  three  miner  leal  solut'oiv  osc'llate 
ahout  the  exact  solution.  The  user  of  my  of  these  approximate  solutions 
would  not  have  anv  automatic  ind’carton  is  to  whether  his  numerical  solution, 
for  any  given  abscissa,  is  above  or  below  the  exact  solution. 

The  above  exanp’es  clearly  demonstrate  the  changing  positional  tola 
tionships  between  approximate  and  exact  solutions  for  even  relatively  simple 
Initial  value  problems.  At  present,  there  are  few  ways  to  determine  the 


nature  of  this  relationship.  One  naive  approach  is  to  construct  a graph 
of  both  the  approximate  and  exact  solutions  for  certain  test  problems. 

In  general,  however,  plotting  is  usually  very  ineffective;  indeed,  accu- 
rate plotting  may  be  difficult  to  perform  il  the  exact  solution  is  very 
complicated.  Frequently,  the  user  employs  a numerical  method  because  the 
analytic  solution  is  inaccessible,  or  he  has  little  knowledge  of  its  exact 
form,  thus  making  plotting  impossible. 

Mo s t approximation  methods  provide  a theoretical  error  bound  for  the 
maximum  dist  mce  the  exact  solution  can  deviate  from  the  numerical  solu- 
tion. I'nf ortunate ly , most  of  these  error  bounds  are  defined  'or  the  abso- 
lute value  of  the  deviation  so  that  the  user  can  sti’l  not  be  certain  as  to 
which  side  the  exact  solution  resides  with  respect  to  the  numerical  solu- 
tion. More  importantly,  these  hounds  are  usually  difficult  to  evaluate  in 
a practical  manner  and  often  require  the  introduction  of  many  simplifying 
assumptions.  Even  when  they  can  he  evaluated,  these  bounds  are  usually  not 
optimal  and/or  do  not  take  into  account  the  floating-point  roundoff  error 
propagation  (which  can  be  substantial  in  some  problems)  encountered  when  an 
algorithm  is  implemented  and  executed  on  a computer. 

From  the  above  discussion,  it  should  be  apparent  that  in  many  situ- 
ations the  user  of  numerical  methods  ordinari.lv  has  no  efficient,  automatic, 
and  practical  indicators  o(  the  positional  relationship  between  approximate 
and  exact  solutions.  One  way  to  establish  an  indicator  is  to  develop  nu- 
merical methods  which  provide  upper  and  lower  bounds  to  the  exact  solution; 
techniques  which  bound  the  exact  solution  from  above  and  below  are  called 
"bilateral"  (two-sided)  methods.  One  such  bounding  approach  is  known  as 
interval  analysis;  this  technique  produces  upper  and  lower  bounds  to  the 
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solution  of  large  collections  of  scalar  and  matrix  mathematical  problems. 

Before  examining  the  properties  o'  interval  analysis,  we  shall  look 
briefly  at  the  engineering  applications  and  other  advantages  offered  bv 
the  use  of  bilateral  methods  in  general.  Even  though  bilateral  algorithms 
have  not  as  vet  been  widely  available,  several  applications  have  already 
been  described  in  the  literature.  • - example,  Baranov  (ll>o-*l  reports  that 
nonlinear  differential  equations  representing  flexural  and  torsional  vi- 
bration of  beams  or  buckling  ■?  rods  often  cannot  he  easily  solved  either 
analytical  lv  or  bv  existing  numerical  methods  in  such  a w v »h.  posi- 

tional relationship  between  the  ivorox  .mate  and.  exact  solut : ons  i.  cu- 
rate’. v revealed.  A bilateral  approach  is  helpfi.’  here;  ar.r-ov  ' deve- 
loped such  a method  for  hie.  specific  oquat  • ons  . In  r.av.i  ’ warfare  -•  • clems 
involving  the  use  of  ’.audios to;  equa’  :ous,  Vobrv  ( 1 ' as  d.o’  i ne,  two 

sided  bounding  scheme  which  offer  a means  of  studying  variations  o’  -er- 
tain  parameters  before  the  equations  are  solved  *>v  sea  convent  Iona  1 .non- 
bounding) method.  Bo ley  (19bo)  has  used  a bilateral  approach  t or  problems 
of  heat  conduction  in  melting  or  sol  difving  slabs;  his  work  car  ilso  pro- 
tpproximatlons  for  certain  types  o!  lerodynatnlc  »lat  >n 
Appl  and  Hung  (19b-*)  have  applied  • two-sided  technique  to  . 'ntvv.  u s ecu  - - 
1 ibrium  nroblems  such  as  a fin  temperature  problem  in  which  there  inter 
nal  beat  generation.  Ispolov  and  Appl  (1971'  hive  employed  i bllarera 
method  for  i problem  of  self-sustained  vibration  o*.  an  uitonomous  system, 
kahan  (l"'.’)  has  developed  an  ellipsoidal  bounding  technique  which  can  he 
applied  to  the  N-bodv  problem.  Pwo-sided.  bounding  methods  have  .t’so  been 
employed  bx  Weinstein  and  Stene.er  (ll>'.')  tot  problems  involving  vibrations 
of  cantilever  plates  or  for  energy  levels  in  quantum  mechanics.  In  most  of 


the  above  examples,  no  general  techniques  have  been  devised,  lust  methods 


to  solve  very  specific  problems. 


In  addition  to  the  aforementioned  engineering  applications,  there 
are  several  other  benefits  to  be  gained  from  the  utilization  of  two-sided 
algorithms  for  the  solution  of  problems  in  numeric  il  analysis.  Kor  example, 
the  bounds  produced  by  a bilateral  method  can  provide  information  ibout  the 
region  in  which  a solution  res:  ; .hi-  max  be  sut  f ic  -.  out  foi  some  appli- 

cations (e.g.  see  Willson  (1968)1,  im'iu.’.iiv.  the  solution  of  algebraic  equa- 
tions, ordinary,  and  partin'  differential  equations,  and  eigenvalue  problems. 
Acquired  knowledge  of  the  solution  region  via  efficient  bilateral  techniques 
could  be  of  assistance  in.  a variation  of  parameter  study:  sue  was  the  moti- 
vation -f  Fabry's  nested  bound  appv  \ eh  (1°  dl  . A so,  the  ins  * eh  t acquired 
about  the  solution  space  could  aid  n selecting  the  most  appropriate  con- 
ventional (nonhound ing)  nuraei ‘cal  method  to  use  on  a set  ot  problems  having 
solutions  in  the  same  or  neighboring  regions.  Furthermore,  good  starting 
values  for  these  numerical  methods  cov’  d be  produced  by  two  -sided  algorithms. 

A bilateral  approach  can  he.-  in  error  estimation.  The  error  of  a 
numerical  solution  obtained  by  conventional  techniques  can  be  compared  with 
the  bounds  obtained  bv  a two-si  c.-'d  scheme,  thus  providing  an  expltc  ^ esti- 
mate of  the  erro>  associated  with  the  approximate  solution.  Also,  the  aver 
age  of  the  upper  and  lower  bounds  produced  by  a bilateral  method  can  serve 
a-,  an  improved  approximation  to  the  exact  schit  on,  in  which  case  this  new 
approx  rat i on  is  in  error  bv  no  more  than  the  absolute  value  of  half  the  dif- 
ference between  the  utilized  upper  and  lower  bounds  (it  it  can  bo  assumed 
that  propagated  co-vputat  icinal  erroi  in  calculating  the  bounds  is  relatively 
small).  Thus,  bilateral  methods  not  only  guarantee  a priori  the  relative 
positions-  of  their  Iterates  with  respect  to  the  exact  solution,  hut  also 
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assist  in  providing  an  approximate  solution  of  known  accuracy.  This  afore- 
mentioned trait  is  quite  significant,  since  most  currently  available  methods 
cannot  offer  an  explicit,  easily  computable,  and  precise  estimate  of  the 
error  (theoretical  or  actual)  associated  with  their  approximations. 

In  our  above  comments  concerning  engineer i„g  applications  of  error 
bounding  *ethods.  we  gave  a few  examp l < of  the  a istance  a bilateral  ap 
proach  offers  when  dealing  with  differential  equations  which  arise  from 
nanv  phys  - i pi  blems;  no  dc  bt  turner  • ,t  et  ■ icatio  , ■ . 

VerV'  A"  trr°r  Sound'nR  approach  should  be  considered  when  he  exact  solu- 
tion is  quite  difficult  .or  Impossible)  to  obtain  or  easily  represent,  and 
uhen  no  -,  standard  numerical  techniques  do  not  clear!,  reveal  the  relative 
position  of  the  approximate  and  exact  solution..  Th . s was  indeed  the  situ- 
: n*  “ rted  ' 1964)  f ranics  problems.  In  - . 

aref  ° " ' " tHe0r>  SCral  ov  ; tn  - the  exact  so 

.ut  on  is  known  only  for  relatively  si-,’0  . • 

■ ^ • tn  comp  .ex  cases,  an 

upper  or  ’over  solution  bound  could  he  very  helpful.  Information  ^Pu!r  the 
solution  region  obtained  from  a bilateral  approach  mav  be  beneficial  n de- 
signing a piece  of  equipment  (e.g.  in  defining  a d-v  ee's  rang,  of  iim  ta- 

tl<M0'  ’ 1ClCn  °TOpUter  ' « >n  . tw  - led  sch«  can  b<  ■ ■ 

m a quantitative  analysis  of  parameters  and  initial  conditions ; such  an 
application  mav  produce  competitive  or  o'gni-canr'v  better  results  than 
th°8e  0bMtned  by  COnw”tional  technique-..  From  the  above  discussion  ir 

Sh°Uld  hP  r°fldiIV  apparent  that  ^ advantage  offered  to  engineers  and 

scientists  is  a possible,  practical  way  of  avt omatical.lv  assessing  the  posi- 
tional relationship  between  approximate  and  exact  solutions.  Thus,  interval 
analysis  is  one  means  of  providing  such  a tool. 


Now  wo  proceed  to  discuss  some  of  the  specific  attributes  of  interval 
analysis.  The  reader  interested  in  extensive  coverage  of  interval  analysis  is 
urged  to  consult  Moore  (1966);  details  relating  to  the  MRC  interval  analysis  pack- 
age are  covered  here  and  in  greater  depth  by  Yohe  (1977)  . We  focus  attention  on 
some  of  the  elementary  properties  of  interval  analysis  and  list  without  proof 
(proofs  can  he  found  in  Greenspan  (1971)  and  Moore  (1966))  some  of  the 
relationships  which  users  of  the  interval  package  should  be  aware  of.  When  in- 
terval analysis  is  applied  to  a mathematical  method  it  generates  bounds  which 
form  a rectangular  parallel? iped  which  encapsulates  the  exact  solution  to  a prob- 
lem. The  technique  utilizes  interval  versions  of  the  primitive  arithmetic  op- 
erations and  combines  these  intervals  using  certain  rules  (defined  below)  to 
decompose  mathematical  expressions  and  problems;  such  decomposition  is  applicable 
to  an  extensive  range  of  numerical  problems. 

The  intervals  we  are  dealing  with  are  closed  and  denoted  by  [a,b]  where 
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[a,b]  ={x  i a<x<b}. 

The  elementary  rules  of  arithmetic  for  intervals  are  as  follows: 

[ a , b ] + [ c , d ] { x+y  | xe  [ a , b ] , ye  [ c , d ] }=  [ a+c  ,b+d  ] 
l«»b]-[c,d]  = {x-y|xe[a,b] , ye[c,d]  }-[a-J,b-c] 

[a,b]  • [c,d];{x'v|xe[a,b],  ye [c ,d  ! 1 = [min(ac ,ad ,bc ,bd) ,max(ac,ad,bc ,bd) ] 


If  Of1  [c,d]  then 


+ *=  f--1  xe  [a,b]  , ve  [c  ,d  ]'■=  rmin  (*}•'•  !•-) 

it,'1  y d c d c c c a c 

From  the  above  rules  It  can  also  be  inferred  that 
- [a,b]  » [-b,-a] 

I - fi,-]  if  Oe'  [c  ,d  ] . 
lc ,d ] d c 

Several  theorems  which  are  analogous  -o  those  for  the  real  line  can  be  establisbe. 
for  intervals: 

commutative  law  of  addition  [a  ,b]-*-'c  ,d]  = [c  ,d]  + [a  ,b] 

commutative  law  of  multiplication  la, hi  • [c,dl=*[c,cU  • [a.bl 

associative  law  of  addition  ([a,b]+[c,d])+[e,f J-[a.b J+([c,dI+[e,f 1) 

associative  law  of  multiplication  ([a,b]  • [c,d]l  • [e,f]=[a,b’  ■ ([c.d]  • [e.f]) 

Ar.v  real  number  a,  his  an  interval  (degenerate'1  counterpart,  namelv 
[a, a].  r,’us  the  interval  [0,0]  corree  onds  to  0 and  f 1 , 1 ’ corresponds  to  1 and 
hence  [ 0 , C 1 serves  as  an  additive  interval  identity  (i.e.  [a,b]+[0,0]=[a,b])  end 
[1,1?  is  a multiplicative  interval  identity  (i.e.  [a.b]  • [1,1] - [a ,b) > . Further- 
more, the  above-mentioned  rules  for  addition,  subtraction,  multiplication,  and 
division  of  intervals  when  applied  to  degenerate  intervals  lead  to  obvious  in- 
terval results  analogous  to  the  non-interval  results: 


[a, a]  + [b,b]  - [a+b,  a+b] 
la, a]  - [b.b]  - [a-b,  a-b] 
[a, a]  • [b.b]  - [ab,ab] 


U.a]  „ f £ a 

[b.b]  lb’bJ 


for  b i 0 . 


Extensions  of  interval  operations  to  functions  are  also  possible.  For  ex- 
ample, Moore  defines  for  f(x^,x2,...  x^)  a real  rational  expression  in  which 
each  independent  variable  x^  occurs  only  once  an  interval  extension 

F(X1,X2,...  Xn)  = {f(xrx2 xn)!xieXi,  i = 1,2, ...n} 

where  each  x^  ranges  independently  over  the  corresponding  interval  . 

Additional  properties  of  interval  analysis  can  differ  somewhat  from  our 
experience  with  non-interval  arithmetic.  For  instance,  the  definition  for 
[a,b]n  for  n a positive  integer  is 


' [an,bn]  if  a > 0 

[a,b]n={xn|xe[a,b] } *=  [0,max(an,bn) ] if  Oe[a,b]  and  n is  even 

[an,bn]  if  Oe[a,b]  and  n is  odd 
[bn,an]  if  b < 0 and  n is  even 

< [ an  ,bn ] if  b < 0 and  n is  odd 

But  [a,b]n  is  not  necessarily  equal  to  n factors  [a,b]  • [a,b]  ....  [a,b] 
e-g.  t-1,1]2  - [0,1]  but 

[-1,1]  • [-1,1]  = [-1,1]. 

Another  example  is  given  by  the  lack  of  >n  exact  counterpart  of  the  distributive 
law  for  real  numbers;  instead  we  have  subdistributivity 
[a,b]  ([c,d]+[e,f  ] )<^-[a,b]  [c,d]+[a,b][e,f] 
i.e.  exact  equality  does  not  necessarily  hold  for  all  intervals.  If,  however, 
for  any  real  number  t or  for  [a,b]  • [c,d]  > 0,  then  the  distributive  law  does  hold: 
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t(l«,n  • u.h])  - ♦ tlg.h] 

[e,f]([a,b]  ♦ (c.d))  - |*,f][a,b]  ♦ [a.fllc.d] 


Also  Interval  arithmetic  exhibits  t be  propatty  of  inclusion  monoton  lc  1 ty , l.e. 
for  [a.b]  ' (c,d)  and  [e,f].  [g.h| 

then  [a.b]  + (c.f)  t [c.d]  *■  [g.h] 

[a.b]  - [e.f]  C.  Ic.d)  - [g,h] 

[a.b]  • I e , f 1 C [c.d]  • [g.h] 

[a.bl/le.fl  cr  [c.dl/ls.h]  If  0/[e,f]  . 

From  the  rules  for  Interval  arithmetic  it  is  apparent  that  the  bounds  for 
the  computed  intervals  usually  involve  several  arithmetic  operations  which  do 
not  necessarily  produce  exact  results,  especially  when  dealing  with  finite 
floating-point  representations  as  used  on  most  machines;  this  error  in  computing 
the  error  bounds  requires  special  attention  to  ensure  that  the  endpoints  of  the 
interval  resultant  are  specified  so  as  not  to  exclude  values  (including  the 
exact  or  true  result)  that  should  bo  included  in  the  interval  bv  the  above- 
mentioned  rules  for  interval  arithmetic.  This  dilemma  is  resolved  by  using 
so-called  directed  rounding  with  all  interval  resultants,  i.e.  the  left  end- 
point  is  rounded  down  and  the  right  endpoint  is  rounded  up;  further  details  are 
given  by  Tohe  (1973,  1*577).  Greenspan  (1971)  provides  a simple  decimal  example 
of  this  situation.  Assume  that  a computer  can  only  store  4 decimal  digits  for 
each  intervat  endpoint  but  that  the  machine  determines  that  the  result  lies  in 
the  interval  J ■ [ -11127.  •.’11231.  If  symmetric  rounding  to  four  decimal  places 
is  used  on  the  endpoints,  then  we  have  3^  « [*1113,  *2112]  hut  we  note  that  a 
possible  true  solution  x contained  In  3 need  not  he  contained  in  J ^ because  3 
is  not  contained  in  Jj.  Ity  using  the  ahove-ment toned  directed  rounding  on  the 
interval  3 we  obtain  J.,  " | •1112,  •2113]  and  note  that  3 j does  indeed  include  all 
values  from  Interval  3.  Of  course,  when  applying  directed  rounding,  one  should  n 
unduly  widen  the  interval  bovond  the  amount  necessary  to  compensate  for  the 
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roundoff  error  propagation  or  otherwise  the  Interval  widths  can  grow  un- 
necessarily large.  Thus  even  though  in  the  above  example  directed  rounding  was 
applied  at  the  decimal  digit  level,  in  practice  directed  rounding  is  applied 
at  the  bit  level  and  is  dependent  on  the  arithmetic  pecularities  of  the  floating- 
point arithmetic  on  the  host  machine  on  which  interval  analysis  is  performed. 

Further  consequences  of  using  interval  arithmetic  are  discussed  in  Section 
5 of  this  paper. 
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APPENDIX  C:  STANDARD  FORTRAN  NUMBER  AND 
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SIGN 

INTEGER 
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EIXEDPOINT 

EXPSEP 

EXPONENT 

NUMBER 


EN’DPTSEP 

COMMA 

INTERVAL 


0|l|2|3|4|s|6|7|ti|9 

i-  ♦ !- 

I-  null|<sicn>|<intecer><digit> 

I"  . 

I-  <INTECERXRADIX>I  <PIXEDPOINT><DIGIT> 
l “ F.  | D 

i-  <SIGN>  1 <EXPSEl"  j <EXPSFP><SIGN>  | <EXPON£NT> <DIGIT> 
l-  < INTEGERS  J <FIX.sr>POINT>  | <INTE  CERXEXPONSNT>  [ 
<FIXEDPOINT>  <EXPONENT> 

I - l 

1“  t 

I-  <N'JMBER>  | (<NL’MBKR>)  f <NUM3KR><ENDPT5EPXNUMBER>  [ 

(<NUMBER><ENDPTSEr><NUM3ER>)  ] ( <NUM3ER» <COMMA> <NUMBER>) 
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